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POREPORD 


Ihls technical memorandum has been prepared as a documentation for an 
objective analysis scheme developed at the Laboratory for Atmospheric Sciences. 
Tlie analysis procedure was developed for numerical prediction studies with 
the POGE data In an assimilation mode with the 4th order 4* x 5* GLAS GCM. 
Although the objective anal"sls scheme is continually evolving this documen- 
tation should be a useful iv^erence. 


layman E. Baker 
1 March 1981 
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1. WroODUCTICU 


Preparation of the Initial conditions for mjmprlcal forecasts has hecomt* 
Incrraslnfcly complex with the augmentation to the conventional database (sur- 
face, rawlnsonck>, pilot balloon, and aircraft reports) of data from spaceborm' 
observing systems (e.g, temperature retrievals from TIR06-N, cloud-notion wlml 
vectors from gmistatlonary satellites''. Many of the non-convent tonal observing 
systenici provide temperature, wind, or pressure only (see Table 1) with dif- 
ferent ermr cliaracterlstlcs. This variety of data and errors must lie ad- 
dressed In an objective analysis scheme. 

This docuna»ntatlon describes an objective analysis procedure developed 
at GIAS for the assimilation of the heterogeneous P3GE database with the 2nd 
and 4th order general circulation models. The objective analysis schemi* Is 
a modification of the Grossman (1059) scliene based on a metlKxl developixl \yy 
Hergth^trsson and nocis (1955). Additional Information on the applications of 
the olijivtlve analysis pnx'rtxlure in numi»rlcal priHllctlon studies with the 
RX'iK data may lx* found In Waken (1981). The gent'ral circulation nxlel ((X.31) 
used most frequently In these studies, the GLAS 4® x 5® 4th order model. Is 
described In detail In Kalnay-Rlvas et al. (1977) and Kalnay-Rlr'as and 
Hoitama (1979a). nocumentatlon of the model Is available In Kalnay-Rlvas 
and Hoi tana (1979b). 

Section 2 outlines the preparation of the POGE data for analysis. The 
ccmponents of the objective analysis scheme are discussed in Section 3. 

Section 4 contains the organl7Atlon of the objective analysis program and a 
d«>scrlptl«x» of the required data sets. Finally, the program logic and a 
detailed description of oaclr subroutine is pt>'st*ntod In Section 5. 


Tablo 1. ntita utMlxod from the P3GE observing systems 



2. PREPARATION OF 'HIE FOGE DATA FOR ANALYSIS 


For this study the official FGGE Il-b database, acquired from Sweden, Is 
utilized. A preprocessor was developed to connect or delete "well-defined" 
errors and to order the data by latitude and longitude. The modifications 
to the database are described below. 

All duplicate reports were eliminated. Fixed reporting stations were 
checked against a list of stations (provided by P. Kallberg of the European 
Centre for Mtxlium Range Weather Forecasts) kncjwn to have had at least one 
occurrence of an Incorrect geographical location. All ship locations were 
checked for possible erroneous land coordinates. It was not possible, however, 
to verify the correct position of ships over the ocean. Errors In excess of 
100 mb in the sea level i.essure reports (ship or land), caused by tiie ambiguity 
In the? report code, are corrected In the objective analysis program by adding 
or subtracting 100 mb. We typically find no more than 5 to 10 such occurrences 
In a 6 h period. 

All ASDAR temperature data were deleted because of an excessively warn 
bias found in the analysis of that data. Rawlnsonde temperature and height 
data were corrected at 100 mb and above for the effects of incoming solar 
radiation on the thermistor with routines provided by G. Costello and J. 

Laver of the National Meteorological Center (NMC). Microwave retrievals 
from TIR08-N with preclpitable water contamination were eliminated according 
to the criteria of Phillips (1980). 

All cloud-motion wind vectors reported at 700 mb or below are re-asslgned 
uniformly to 900 mb In the objective analysis program. We also re-asslgn 
all cloud-tracked winds reported above an estimated model tropopa.use to the 
tropr^Ause level. Cloud-motion winds reported between 400 mb and the estimated 
model tropopause are re-asslgned to the pressure level whose model temperature 
corresponds to the reported cloud-top temperature. 
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^ M O&JBCT WE A NALYSIS SCHEME 

In this section wo describe the current version of the objective analysts 
and assimilation procedure, Hio objective analysis scheme is a successive 
corrections method (SOU) of the Cre^nan (1950) type which involves successive 
modifications to the first guess fields provided by the 4th order GI^AS cxrw 
(Kalnay-Rivas and (ioitsma, 1979). Bustward and northward velocity corpononts 
u and V, geopotential height z, and relative humidity RH are analyzed on 
mandatory pressure levels. Surface pressure and temperature are reduced to 
sea level and analjv^ there. 

The analysis of geoi»tentlal height rather than temperature has the 
following advantages: 

1) Significant level terperature data are utilized in computing rawinsonde 
mandatory level heights. 

2) Chlculatton of the gocjstrcphlc wind is straightforward with a helgtit 
analysis, 

:i) Ihe mean orrtjr for large thicknesses computed fix»n TIHOS-N’ infrar\^1 
retrievals is quite small. 

Ttie assimilation procedure involves the Intermittent analysis of batches 
of data grouped in a +3 h window about each sjmoptic time. No attempt is made 
to balance the mass and motion fields. 

3.1 Vertical interpolation between o and P 

The model first guess fields are updated every 6 h at the model o levels 
wliile the data are analyzed at the mndatoiy pressure levels. This requires 
a double interpolation between o and p. In regions of ’■ ’■'f nxiuont updating 
this results in excessive sioothing, ;is McPherson f^t a/U (1979; point out. 

We have attaipted to reduce the effects of the vertical interpolation by 
interpolating only tfie difference between the moilel first c-’css and ‘h'' 
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analyz/ed fields (e.g. Au, Av, ARH) leather than tne analyzed fields 
themselves Oensequontly , if a particular gi'ldpolnt Is not affected by data, 
the dlffti-ence between the first guess and the analyzed field is zero, and 
the first guess is returned to the model. Interpolating the difference 
between the model first guess and the analyzed field also helps to reduce 
Interpolation errors where updating lias occurred. However, even If the 
analysts and updating were performed In the model o coordinate, the 
vertical Interpolation problem would not be eliminated as \fcPherson et al. 
note. 

The following vertical Interpolation procedure is used. Fbr the 
wind analysis, the first guess u and v wind components are interpolated 
(linear in log p) from the model a levels to the mandatory p levels. After 
the analysis is completed, the differences between the model first guess and 
the analyzed fields are obtained at the a levels by linear in log p Inter- 
|X)latlon. To avoid extrapolating above the top model o level (nominally 
Of) mb) to obtain a 50 mb first guess wind field, a 50 mb wind analysis is 
not perfonrcd. Instead, we calculate 

A 

^65 “ * ^'^70 (1) 

V65 ** ^65 ^^’70 (2) 

where U05 and V05 are the model first guess u and v wind components at the 
top 0 level and AU70 Avyo represent the dlffere've in '.’v> mcxlel first 
guess and the analysis at 70 mb. 

The vertical interpolation between o and r, for the '’'^•lat'' v-' humidity 
analysis Ls performed similarly exc;ept that ndative humidity Is first 
calculated from the model specific humidity q at the ^ level'- using thi' 
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approximation 


RH » q(p-0.378es)/(0.622%) 


(3) 


where % is the saturation vapor pressure. Vte then evaluate eg using the 
empirical formula of Tetens (1930) as nedified by Murray (1967) 

Os ■ 6.11 exp [a(T-273.16)/(T-b)l (4) 

where a *» 17.269 and b » 35.86 if T > 273.16 and a « 21.874 and b « 7.66 if 
T < 273.16. Next, RH is interpolated (linear in log p) to the p levels and 
analyzed at 300 mb and below (see Fig. 1). RH is then interpolated to the 
model 0 levels in the same manner as was done for u and v, followed by the 
calculation of q using Bq. (3). No change is made in the first guess specific 
humidity field above 300 mb. 

We analyze the geopotential height z of the mandatory p levels, whereas 
a first guess temperature field is provided at the model a levels. Tliis requires 
a somewhat more complicated procedure tlmn was needed for the wind and moisture 
analysis. To prepare the model first guess z field for updating with the 
observations at the mandatory p levels, a linear in log p interpolation is 
performed to obtain the model o-level teimjeratupes^ at the mid-mandatory p 
levels. Hie hydrostatic equation 

p-l3p/te =» -g(RT)-l (5) 

is then integrated from sea level to 50 mb to obtain the model first guess 
geopotential height at each mandatory pressure level. In Bq. (5), g is the 
gravitational constant and R the ^s constant for dry air. 
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Virtual temperature is actually interpolated. 
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First, to compute the first gu^s geopotential height of the 1000 mb 


surface z we Integrate Bq. (5) from sea level to 1000 mb to obtain 

1000 * 

^ 1000 “ /lg/R(*n(Pst/1000)) + 0.5B] (6) 

wliere 8 = 6.5 °C/km, the standard atmospheric lapse rate, Tg£ the sea level 
t{^nperatui"e , and Pg£ the sea level pressure. In Rq. (6), Tg£ is calculated 
as 

with Zg and Tg tlie elevation and temperature, respectively, of the model lower 
boundary (0 * 1). The sea level pressure Pgj is determined from 

Pan Pgexp{gZg(RT^)"^ ] (8) 

\^ 11 cre pg is the model surface pressure and Tm the layer mean temperature 
lx-* tween sea level and the surface expressed as 


T,n = Tg + 0.5teg 


(9) 


To calculate the first guess geopc^tentlal height at the n levels a’x>ve 
1000 mb, we integrate Bq* (b) from 1000 mb to 50 mb using 


^^k+1 = ''ic (Pk - Pk+l)/Pn 
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where k indexes the mandatory levels and l the mid-mandatory levels, Cp Is tha 
specific heat at constant pressure, and < » R/Cp. Tlie raid-mandatory level 
pressure p^ raised to the power < in the denominator of the bracketed quantity 
in Bq. (10) is computed as, 

pj - <pfi-P{:})/(Pk - Pk*i) ai) 

which assumes a constant potential temperature in each layer after Phillips 
(1974), p£ is defined similarly in the forecast nodel (Kalnay-Rlvas et al. , 

1977). In the interpolation of the o-level temperature profile to obtain 
in Rq. (10), described previously, the pressure pjj corresponding to each Tii 
is computed using Bq. (11). 

After the geopotential height analysis has been completed, mid-raandatoiy 
level tanperatures are calculated. Using Bq. (10) we solve for such that 

T’ji = BCp'^(^+l-^k> (Pk - Pk+l)"^ (12) 

The difference between the analyzed and the first guess is then interpolated 
linear in log p to the o levels and culded to the original model temperature 
profile. The updated o-level temperature To may be expressed as 
A 

To = To + ATo (13) 

A 

where To is the first guess temperature and ATo interpolated difference 
between the analyzed T and the first guess T. We assume that the difference 
in dry bulb temperature ATo veiy nearly equal to the difference in virtual 
temperature. 


3.2 Data selection and horiaontal consistency 


After the observational data liave been Ingested, a comparison is nude 
between the data and the model first guess fields interpolated to constant 
pressure surfaces (see Fig. 1) as previously described. Data which differ from 
the model first guess by a specified amount are flagged as suspect. Hie im.ximum 
allowable difference between the observations and the model first guess are 
shown in Table 2. The choice of these "acceptance criteria" is arbitrary. 

We were guided in our selection of these criteria by the desire to use all 
available "good" data at the risk of including a few "bad" reports. 

After all the data have been checked against the model first guess, those 
dotn flagged as suspect are rechecked by comparing the average difference between 
the nearby reports (observations within a 5* I’adius) and the model first guess 
w^ca the difference between the model first guess and the observation in question. 
An observation is rejected whose difference from the model first guess differs 
from the average of tlie nearby differences by more tlian the limits given in 
Table 2. 

3.3 Interpolation of observationaJ. data to the tiorizontal grid 

Hie method of interpolating the observations to the 4® x 5® analysis grid is 
based on the Cressman (1959) scheme, but modified to treat variable data density 
and quality and the first guess accuracy of tiie prediction model. Tlie interpo- 
lation procedure is described below. 

Let the total correction to tlie background field be expressed as 

N M N 

Ck = )' WiAFt/( y/i * I Wi) (14) 

i=l j=l i=l 


to 








wliert> k Indexes the number of scans (k -1-3 In the present scheme), I the N 
obsfirvatlons available In the kth radius about a model grldpolnt, J Uie 
surjxKindlng M grldpolnts about the grldpolnt that Is to be updated, and 
denotes the deviation of each observation from tlie first guess field. Tiie 
deviations art* calculated using a bilinear interpolation of tiie four model 
grldpolnts surrounding each observation. Tiie iiackground field Is tlien subtracted 
from the observed data. 

In Fxj. (14), and wj represent the weight given to tiie Ith c)hsorvatlon 
and Jth grldpolnt pjvdictlon, respxjtlvely. We define 

'*'1 = 


and 


Wj “ qj(Kg - r-j)/(R^ rj) (U*) 

wlH*re R|^ Is tia* size of the kth radius of Influence abait eficii jnodel grldgolnt, 
Uie dlstfutce from the ith observation to ti»e model grldpolnt to be update!, 
juid rj the distance frar the jth surnxmdlng grldpolnt to the grldpolnt 
Ixjlng analyzed. %. determined as a function of data density. Is defined as 

Rk = Okd (17) 

where d Is tiie avera^;e distance between observations within a glveti radius 
iibout eacii grldpolnt. In Hq. (17), d nvay be exivx*sseci !is 


d = (!i/N)^/2 


(18) 


where N is the number of (Nervations within a radius r^j (800 km in this study). 

In Bi. (17), tlie coefficient cj^ is chosen to be 1.6, 1.4, or 1.2 for scans 
1,2 or 3, respectively, after Stephens and Stitt (1970). 

In Bqs. (15) and (16), and qj are pre-assigned obsen'ational and 
first quess weights, respectively. Values for represent the ratio of the 
mean ^uared error of the rawinsonde data to the mean squared error of all 
other data. Rawinsonde data are given full weight in this analysis scheme 
with respect to observational accuracy. For non-rawinsonde data qj^ mnges 
from 0.0 to 1.0. In cases where the mean squared rawinsonde error is greater 
than that for other data, q^ ** 1.0 for tliat data. 

Tlie nns differences that are used to cxmpute the observational weights 
are (xmputed using the NMC final analysis. A +^3 h window is used about each 6 
h analysis with the same data checking criteria given in Table 2. Tie matchups 
between the NMC analysts and the observations are acccmplished with a horizontal 
bilinear interpolation of the NMC analysts to the l(X»tlon of each observation 
at each mandatory level. For ofi'~level wind reports, the NMC analysts is 
interpolated verttcjally linear in log p to the reported pressure levels. 
Differences between the analysts and the observations are then obtained and 
assigned to the nearest mendatory level for inclusion in the rms cmlculatton. 

The ms differences for off-level aircraft temperature data are obtained similarly 
except that geopotential height is computed at the nearest mandatory level for 
the aircraft report. 

The quality of tlie first guess field qj represents the ratio of the mean 
squared error of the rawinsonde data to the mean squared error of the 6 h 
prediction augmented by the update history of each gridpoint at each level. 

Tie mean siiuaixjd prediction error is not pemitted to gi*ow beyond 72 h in 
regions of infrequent updating. No attempt is made to account for varying 
data quality in the update history (NtePherson et al . , 1979). In compiling th** 


error statistics for the model prediction, each G h analysis during a 32 d 
assimilation experiment utilizing the PQGE Level Il-b data Is used to verify 
the 6 h first guess for 127 cases. 

For the horizontal interpolation of the wind data to the analysis grid, 
the distance weighting in Bq. (15) is modified to allow upwind or downwind 
observations to have greater welglit than crosswind observations os was done by 
Fiidllch and Mancuso (1908) for analyses of Jet maxima. In the present schane 
we mcxilfy tiq. (15) such that 

Wj^ = q^[R^os<> (R^ - r?)sin>{)]/(R^ + r^) (19) 

wliich Is similar to the weigliting function used by Bergman and Carlson (1975) 
for analyses of tropic-al cyclones. 

3 . 4 l^atln g th e Model Fields 

3.4.1 Analysis of the mass f ield 

After first guess fields of Tg£ and Pg£ iiave been obtained using Eqs. 

(7) and (8), respectively, the backgrcxind fields are modified by the oljservational 
data u.*- ing Bq. (14). All available land and sliip reports are analyzed. In the 
present scheme, wind observations do not affect the sea level pressure analysis. 

An 8th order Sl\apiro (1970) filter is applied once after the third scan In the 
analysis of pg^^ and Tge* 

Hie model surface pressure is updated by solving for pg in Eq. (8). The 
vertical coordinate of the model is then adjusted using 

0 =(p - P )/(Ps - P ) (20) 

T T 

v4iere p = 10 mb, the pressui*e at the top of the mcriel. 


3.4.2 Updating the thennal field 


In order to update the mxlel o-level ten^ratureti. a bachf round height 
field Is obtained at the mandatory pressure levels as previously described. 
Height profiles are constricted for singlo-level aircraft tempemture data 
using I3q. (10). A first guess height profile is assumed below the level of 
the aircraft report, with the aircraft temperature data modifying the helglit 
profile at all mandatory levels above the aircraft level. 

Hydrostatic consistency is maintained In the vertical by chcjcklng the 
static stability of each layer against the limits given in Table 3, where D 
represents the actual height minus the reference height and S the layer mean 
stability. The superscript 6 refers to the 6 h model first guess. The 
infonnation for the first 10 levels contained in Table 3 was provided by 
J. McDonell of the NMC Automation Division, Two additional levels were added 
at 70 mb and 50 mb to acccmnodate the analysis at 12 levels. 

First, Zj^qqqIs calculated using Eq. ( 6 ). The SCJi is then applied 
to the 300 mb first guess height field. We followed the procedure established 
in 1971 at NMC (see Tech. Procedures Bulletin No. 63) of selecting the 1000 mb 
and 300 mb surfaces as key levels because of the abundance of sea level and 
aircraft data affecting those levels, respectively. The information at 300 mb 
is reflected in the levels above and below 300 mb through the calculation of 
the first guess (defined in Table 3) at those levels. Similarly, the first 
guess for levels above 1000 mb is influenced by the 1000 mb analysis. 

As an exanple, consider the calculation of D 500 (Z500 ■" reference height), 
defined in Table 3 as 

^500 “ ^^1000 ^300 ^ ^500 
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First, rewrite Eq. (5) as 


3z/3w - -Cp g“^0 (22) 

wtiere i» is the Exner function and 0 the potential temperature. Ekj. (22) may 
tlien be Integrated from 1000 mb to 5(X) mb to obtain 

^500 " *1000 " ®?88o (”1000 ” ”500> (23) 

Integrating Bq. (22) from 500 mb to 300 wb yields 

*300 " *500 ■ ®5o8 ("500 “ "300> (24) 

In E!qs. (23) and (24). represent t!»e mean potential 

temperature of the Indicated layers. We may then use 

f^SOO " ^2500 “ ^72 m (25) 

and Elqs. (23) and (24) to obtain Bq. (21). In Bq. (21), Sqqq is computed 
from the model first guess with the following expression 

^500 " ^^00 " ^000 “ *^ioo 

After calculating D500 in Bq. (21), the SCM is applied to Z500 ^ 
stability factor S500 is computed as 

S5OO " (^00 - (^000 - *617 D300 (27) 

S500 is reset to -204 m or to 125 m for those temperature profiles wtiich are 
either too stable or too unstable, respectively , in the 1000 mb to 300 mb layer 
Dgoo is then recalculated after substituting Sqqq from Bq. (27) for !^qq in Bq. 
The (ialculation of D500 snd S500 Is repeated for t)ie second and third scans 
with t)»e stability checked and reset, if necessary. As with the sea level pres 
sure and temperature, an 8th order Shapiro filter is applied once to the height 
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( 21 ) 


field after the third scan. Iho other levels are analyzed in the order given 
in Table 3. The model o-lovol tcmperatiircs are then updatcid in the nanner 
previously described. 

3.4.3 Analysis of the wind and nplsturo fields 

An analysis of the wind field is porfonned at 11 mandatory levels (1000 mb 
to 70 mb) with corrections made to the backgreamd field using Bq. (14). 'Pio 
30 mb analysis Is excluded to avoid extrapolation, as discussed previously. 
Cloud-motion wind vectors are re-assigned prior to modifying the first guess 
field ns described in Section 2. Single-level wind data (i.e. aircraft, cloud- 
motion winds, etc.) are allowed to affo^t the analysis at the two mandatory levels 
adjacent to the wind report. We use wind data from ships only in the 1000 mb 
analysis. 

Before applying the SCM, a local geostrophic com»ctlon is applied to the 
first guess wind field using a technique similar to that proposixl by Klstlor 
and McPherson (1975). However, the geostrophic correction is computed from the 
change (analysis-first quess difference) to the mass field (surface pressure) 
only. This approach has been tested successfully by Stone et a l. (1973) and 
is in use operationally at the Australian Numerical Meteorological Research 
Center (K. Puri, personal camunication). First guess and analyzed values of 
^1000 obtained using Bq. (6) from which corresponding fields of the 
geostrophic wind are computed. The first guess wind field is then adjusted 
uslt^ 

V = w + [ gf ^ (TZ^qqq - ^1000^^ ^ 

which is Bq. (1) in Kistler and McPherson. Tn Bq. (28), w is the first uuess 
vector wind field and f the coriolis parameter. The adjustment to ^irsl 
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guess wind field is scaled vertically and latitudinally usii% Gqs. (3.10) and 
(3.20) in Bergman (1979). An 8th order Shapiro filter is applied at each of 
the 11 analysis levels after the third scan of the SOA, The model o- level 
wind field is updated as previously described. 

Relative humidity is analyzed usii% Bq. (14) at six mandator)’ levels 
(1000 mb to 300 mb). There is no interaction between the moisture field and 
the other fields during the analysis. As with the other fields, an 8th order 
Shapiro filter is applied once after the final scan in the RH analysis. 
Updated o-level values of specific humidity are obtained using the procedure 
discussed pi’eviously. 


4. PKaiUAi\l UiLlANIZATION AND DATA Sbl' DiidOtlFnai 


^ Objoctive Analysis Modules 

Tliu l‘X)K'fKAN prugruin wtucti i.* *• »ri‘ornk» the objective analysis consists oi 
subroutines contained in throe source modules: 


• GiVtJ.MX.POKTH.ANALYSIS.hXJrtT 

• IDI2HCDA 

• It>I2U2D 


Objective analysis 
Level H-a unpacking 
Level Il-b unpacking 


The object code equivalent of the objective analys. • is contained in 
;ji;o module; 

6 I726bA0l 


4 . 2 Source modules for the 4t ti Order A^el 

llxecution of tlie objective analysis requires on-line interaction witi» 
the 4Ui order model, nie iOKTKAN programs that constitute tue model are 
contained iti five source modules. 

• I?(JIDHI‘’2 Old master. Main and hydrodynajiucs 

• ITUlDHiS Old master; Physics 

• I7GIHADN Old master: Long-wave radiation 

• I761DHTI Old master: Stiort-wave radiation 

• GWS.MI.hXJUT.UTILITY.PDd (itKl*'!’!) Old Master: tbst-Fourier transform 

The current version of the model is obtained via updati'.s to the old 
masters and is contained in object form in two modules: 

*• 17tj24TH4 Object code: Main, liydrodynamics.pliysics.atjd iong-ui)d short- 

wave radiation 

• I7G2HFbT Object code: Fast- Fourier transion.T 
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In addition, a library of routines Is ccmtalned In: 


• I527JLIB 


4.3 S ource Module for the Objective Analysis 
The imdule 

• GWS.MI. FORTH. ANALYSIS. PORT 

contains the following 85 subroutines and 4 entries (Indented); 


ADPINS 

FILLIN 

INUPAR 

RELDEW 

ALTER2 

FIND 

KINDEX 

REIJIUM 

ATEOTA 

FINDHM 

DINDEX 

RINDXl 

ATESTH 

FINDLV 

LEQl 

RINDX2 

ATESTP 

FVARIO 

LOG 

SATa«) 

ATESTQ 

GCDIR 

LOCATE 

SATURN 

ATESTT 

GCDim’ 

LQEU 

SBLIZE 

ATESTU 

GEIOSAD 

LQUS 

SCALE 

ATEsrrw 

GBOSM 

MAPP 

SHAPFL 

ATE?T3 

GBOTOT 

NE»ITPE 

SHUM 

CDGRID 

GETAHT 

NEXT 

SIGTOP 

CHKADP 

GETTEMP 

NMCHIT 

SM0C5C 

CXJOT' 

GETFEGW 

NMCSCM 

SMOOTH 

CUTADP 

GETTWND 

PAIRZ 

SMPSL2 

DFDX4 

GRIDCD 

PARA 

SM2D 

DFDY4 

GTOPOG 

PCAL 

SPEHUM 

DIFFRS 

GTTOT 

POTEMP 

SSCAN 

DUMMYZ 

HITE3«P 

PRESIG 

TAP24 

INSMIT 

HUMID 

PSIOtA 

TERP 

TTHC9A 

INCMN 

PSURFE 

TSFACT 


INSEAL 

PTOSIG 

WINDPR 

INSURF 

The module 
• I612HC9A 

contains the following 9 subroutines: 

PUTWS 

QE 

READ IN 

WIND2 

$SM009 

TTHC9A 

TTINT 

FIX 


GCTCAT 

UPKID 

UPKST3 


LOAD 

UPKSTA 

UPKST5 



The module 


• I612U2B 


contains the following 9 subroutines: 


iXJHM2l4 

UIJK2BW 

IFUXtolA 

UPK288 

KEEU37 

KFWJMA 

Ul^UU 

NEWIiDU 

TBSDAT 


‘4 . 4 Job Contro l La nt^uago l‘or Execut ion 


/ID 


/DD l> 

mNAME=Iiil2UaM 

/DU U 

DUMMY(IJJ) 

/EXl-JC 

bXJKTX 

/DD O 

D8NA.vlE»I(J12U28 

/UD (3 

DUMMY (183) 

/b’XLL' 

FOK'rX 

/JOBLIri 

D8NA*\lE=lb27JLlB 

/Sl/iE='15(X)K 

/DD 8 

D8NAME=»17U24TH4 

/DD 9 

D3NAME= 17U2HFb’T 

/DD 19 

D8NAMl>=I7U2FA0i 

/liXhL’ 

LDK, 'BEX' 


INCLUDE 8 
INCLUDE D 
lt4Q.UDE ID 


/IX) 

D 

DUMMY (2G00) 

/DD 

8 

UEVT-=2409 , HECFiM=U , LKELa^7U80 

/IX) 

11 

DEVIH.XXiE , LRECLr»84 , 8PACE=10 

/DU 

12 

DEV'l'=l\42bO , ltlJCFiNt=U , LRECly=20000 

/UD 

17 

D8I>IAME=lb27SMFZ 

/iX) 

4U 

D8N/lvlE=ib27ALBE 

/DD 

4D 

DSNAMD=CIAPl’7‘JT 

/i)D 

Oi 

D8N/VVU'J-^;i)BY8b’'ri 

/DU 

x> 

D8NA.vU')=l7bLDAl\) 

/DD 

IXJ 

U8NA.MbMJlBKDATl 

/IX) 

21 

D8I4AMt)=I7GaMb'17 

/DD 

23 

D8NA./LE= I7G3AAGG , HbX,>’M=F , LttECL=372 , U80lX4=DA 

/IX) 

24 

DUMMY (2GD0) 

/DD 

2b 

LEVi'=24D9 

/DD 

•2U 

D8NAME=17G3HM8A 

/DD 

Gl 

D8NAMl-.= 17G3t‘V*^3 

/DD 

81 

l)8NAi\«E= I7G3ShM , KECFM=F , BLK8 i/=29G0 , LHECL-=2DGU 

/IX) 

82 

DSNA»vlE=I7G38EAZ , RECFI^F , BLK8 IZ=29GU , UU-;CLr=29G0 

/DD 


U8NAME= I7o3/\/V3B , HECb’ivNF , BJJv81Z=29G0 , LliECIj=29GU 

/UD 

84 

D8NAME=17G3SEA1 , RECF'M=F , BLK8IZ=29GU , UIECl/=29GD 

/DD 

8b 

D8NAMit;= I7G38EAY , RECFM=F , BLK8 1Z=2UGU , L:UXlr=29G0 

/UD 

8G 

D8NAME= I7G3AAGD. HhCFM=F . BLK81Z=29G0 , LRECL=29G0 

/UD 

87 

U8iMA.>ll'>17G3PFr2 

/L-m; 




INUUTZ 

Jil-M) 

JsiAlPUT/, 

JiCNI) 

r 


4 . i) In put D ata Pa rainet ers 
4.5.1 Name list INPUTZ 

Namelist INPUTZ detinus parameters which vary according to the pjirticular 
model experir.ient. Tliey will, therefore, not be described in this report. 


4.5.Z Namelist ALPUTZ 

bX)»{TKAN 

Description 

Default 

'iVp« 

pariuneter 



INSADP 

Time interval for data analysis 

6 

14 

NKOJN’r 

Log 24 tape flag (0 for no output) 

0 

14 

ASTAKT 

Initial conditions flag (1 for o-ievel tape) 

0. 

U4 

MAXDl:X 

Numlier of data types used in analysis 

7 

14 


See Section 4. 5. 2.1 

iOO^.KALSb;. 

\A 

QililAP 

liiiiptical weighting for wind data 

.KALSG. 

Ul 

NSAT 

Number of satellite sounding data sets 

0 

14 

LUdAT 

Ijogical units for sounding data sets 1: 

;10*0, 2;i0*i 

14 


Geostrophic wind correction 

.False 

\A 

WMLAT 

Maxenum latitude for wind analysis 

90. 

114 

NTOPUj 

DSKN for topograpliy used in PTOSIG 

0 

14 


4.5 .2.1 n&ta Type in the (g( array 


location in QX array 


Description 


1 

2 


3 


4-7 


8 


9-11 


12 

13 

14 


15-32 

33-50 


51-56 


Surface, rawinsonde,and pilot balloc« data 

Aircraft 

ASDAR 

NESS EAST, NESS West, EXuxjpean, and 
Japanese cloud-tracked winds 

Drifting buoys 

Wisconsin East, West, and Indian Ocean 
cloud-tracked winds 

Constant level balloons 

NAVAIDS 

Dropwi ndsondes 
TIROS-N soundings 
Special effort soundings 
i'TPR soundings 
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4.U 

8ot iXascription 



[yam 

Status 

Organization 

Description 

8 

Output 

Suq. 

Model nistory ape 

u 

Input/Output 

Seq. 

'IVanporary nuiikjiist cuiv* space 

12 

Input 

Suq. 

0 -level initial conditions 

17 

Input 

Seq. 

Surface geopotontial Ueiglit 

40 

Input 

Suq. 

Surface allx>do 

4U 

Input 

3eq. 

'fOpograpiiy 

bi 

Input 

Suq. 

SST/snow cover data 

50-5(3 

Input 

Suq. 

dadiation data 

21 

Input 

Seq. 

Mandatory- level initial 
conditions 

23 

Input/Output 

D/A 

Temporary storage for 
mndatory-levei daUi 

24 

Output 

Sui}. 

Analyzed iielils at 
liundatory levels 

25 

Input 

Stxi. 

Level I I-H data 

20 

(Jutput 

Seq. 

Values of mis fit 

01 

Input 

Seq. 

(Quality weiglits 

81 

Input 

Seq. 

riuoy data 

82 

Input 

Seq. 

Wisconsin Last and West 
LTW data 

83 

Input 

Seq. 

Wisconsin Indian Ocean 
CrW data 

84 

Input 

Seq. 

Constant level Dal loon 
data 

85 

Input 

Seq. 

NAVA IDS data 

80 

Input 

Seq. 

Dropwindsonde data 

87 

Input 

Seq. 

TIHOS-N data 
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4 .G . I 8. >k.5iul .^istoQr_ 

The inodui tustx>ry tape (16& t’unnaL) is descrluud lu douiil iti 
lilolinann (1980). 

'i**^*^ 9tiHN 1 2 : 8ig>na- LBVsi Initial Oonditions 

Tlio Initial conditioivs for the mcxlul oxpernnentb art* described in 
idelmann (1980). 

4.G.3 U dKN 1 7: S u rface Ueopotential Heitdit 

A description of the surface geopotential height is given in ilaKer 
e t al . (1980). 

4.G.4 OSl^ 40. Ibpogi’aptty 

A description of the topograpijy data.‘set is given in BaKer 
et.al (1980). 

4.G.t3 'Ibmp erature and S now Cbver 


I Uan 
number 

Data eiefnent 

IVpe 

Nuinbt.;r oi 
words 

1 

Simulated tene (TAU) 

K*4 

1 

•} 

Connent array 

K*4 

300 


Snow cover (=1) 

UM 

3313 

4 

Sea surface tduperature (K) 

K*4 

3313 


4.6.6 D8RN 21: Mandatory Level Initial OonditlCTia 


Item 

number 

Data element 

Units 

Dimension 

Type 

l*^ber of 
words 


Record No. 1 





1 

TAU 

h 


R*4 

1 

2 

Spares 




14 

3 

glares 




20 

4 

Sea level pressure 

mb 

(46,72) 

R*4 

3312 

5 

Surface temperature 

K 

(46.72) 

R*4 

3312 

6 

Tropopause pressure 

mb 

(46,72) 

R*4 

3312 

7 

Relative humidity 

(%) 

(46,72,6) 

R*4 

19872 


Record No. 2 





1 

U wind component 

m/s 

(46,72,12) 

R*4 

39744 


Record No. 3 





1 

V wind component 

m/s 

(46,72,12) 

R*4 

39744 


Record No. 4 





1 

Tenperature 

m/s 

(46,72,12) 

R*4 

39744 


Record No. 5 





1 

Geopotential height 

(m) 

(46,72,12) 

R*4 

39744 


4*6.7 DBRN 81, 82, 83, 84, 85 and 86; Level Il-b Efata 


Ttio format for the InternatlcHUil exchange of FtTiE Level II data may be seen 
In the report prepared by tlie Working Group on the Global. Data Processing 
Syston (1979). 

4.6.8 OGRN 25: lovel Il-a Data 


Item 

number 

Data element 

Type 

Number of 
words 

1 

TAU 

R*4 

1 

2 

Year, Month, Day, Hour (GMT) 

1*4 

4 


Pile name (8 elements) 

L*1 

2 

1 

Number of reports 

1*4 

1 


Number of elements in data array 

1*4 

1 

6 

Data array (1) 

1*2 

(Item nunA»r 5) 

A description of the data unpacking routines may 

be seen 

in Edelmann (1979). 


4.6.9 DSRN 40: Surface Albedo 


Item 

Data element 

Type 

Number of 

number 



words 

1 

Sui’face albedo 

1*2 

3312 
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4 .0 . lU UiHN 8 7: aatollite 'fanttwra ture aoun dititi Uata, (NIJSS (Axira tiotmi, 


TIHOS- N and VTPH ) 


Itum 

nuinder 

Data eliaiunt 

'IVpe 

Number of 
woixls 

1 

TAU 

H*4 

1 

2 

Satellite 10 1-TIHOS-N 
2-VTVK 

1*4 

i 

3 

latitude 

rt*4 

1 

4 

Longitude (pcjsltive wtst of Greenwicu) 

K*4 

1 

6 

idem til angle 

K*4 

1 

6 

Elevation of surface (10“^ in) 

K*4 

1 

7 

Skin temperature (K) 

K*4 

1 

a 

Estimated surface pressure (inU) 

K*4 

i 

i) 

Instrument/cliannel combination* 

1*4 

L 

10 

Hetrieval mettiod** 

1*4 

L 

11 

Scares 


3 

12 

'fropopause pressure (mb) 

H*4 

1 

L3 

fropopautse toinperature (K) 

H*4 

i 

14 

Spares 


5 

15 

Ory bulb tenperature (K) 

H*4 

12 

10 

Geopotential height (m) 

H*4 

12 

17 

Virtual temperature (K) 

ix*4 

12 


*icx; » ioooo*v + 1000 *w + ioo*x + io*y + z 

= iOO*X + iO*Y + L 

29 


‘i . o . U X X : j -XiUE i ip ocial fcit'i’ort TlitUj-N Ht^trievai M 


Itura Uuta oiotnuot Type Nunjber of 

nuintxit* Words 


1 TAU 

2 Yeav, Month » Day, Hour ((iWf) 

3 liata soui’ce index 

4 Itopopt identification 

:> Latitude 

() Longitude (0° to 300'^' 

west of Greenwich) 

7 itour of oLsurvation 

8 Minute of observation ( IU“^ ti) 

:) Instrument type 

10 Gpeciai eliort quality ilag 

U Surface elevation (ri) 

13 Surface pressun; imlJ) 

13 Surface ternperature (°C) 

14 Surface dew-point depi-ession ('’C) 

Lo Sea level pressure (mb) 

LG SXin temperature (“C) 

17 Solar ^ietjitii angle -UO® (night) 

to +bO°(day) 

18 Channel cunbination 

ID Standard deviation for iow- level 

ctfantK*l 7 

30 Standttixl deviation for inedimn-levei 

ci;antiei 

21 Value of N'SfAK (.V*’‘iU00) 


R*4 1 

1*4 4 

1*4 1 

R*4 2 

R*4 I 

1*4 1 

1*4 I 

1*4 L 

1*4 i 

1*4 1 

H*4 L 

H*4 I 

K*4 I 

K*4 1 

H*4 I 

H*i 1 

1*4 1 

1*4 i 

1*4 1 

1*4 I 


4.(i.lL OontinuiJd 


Ittjn iJata oLutaunt 'IVp^* Number oi 

tiumbtir woi\ls 


122 

liL>triuvui MettKXl 

1*4 

L 

23 

Orbit liiio number 

1*4 

1 

24 

Joan eiumunt aumbur 

1*4 

1 

25 

Pressure of averatse cloud top (10 nb) 

1*4 

1 

2G 

Cloud cover (%) 

• 

1*4 

L 

27 

Total oi«!one (Dobson Units) 

1*4 

1 

28 

rropopause pressure (10 mb) 

1*4 

L 

2U 

Tropopause temperature (10 “C) 

1*4 

1 

3U 

v^lity indicator for tropopause pressure 

1*4 

1 

31 

Layer mean tenperature ("C) 

H*4 


32 

Layer precipitauie water (ran) 

il*4 

3 

33 

Quality indicator for temperature 
(not used) 

1*4 

l;> 

34 

vjuality indicator for precipitation 
(not used) 

1*4 

3 

35 

Mandatory level ben;Uts (1000-10 mb in in) 

U*4 

15 

3G 

Mandatory level temperature (850-10 mb in °C) 

K*4 

14 

37 

Mandatory level dew point (850-10 mb in ®C) 

U*4 

o 


31 


t>iviisnun Ujt^xc m>d Su bixjuiit io iX^cripuoa 


i).i aubrouUnu Qtilxuti^ ou4uo uco 




^MiA 
TWUITi: 
INC\!N 
h UUAUIN 


-DiNLHiX 

-LHiKlUi 

FijlGIXJP 




-KIND 

"rt/rfcaip 

-SUL WE 

-iHTEfli> 

-KIND 

-ilEUiUM 

-EINUHM 

- KIND 

-iV’AL 

- 1- IND 


— ^Locate 

0JN6Tr 

^^’ATE 


£^'ATE 


— j~ U UCATE 


GEOSAL) 


WArtiO 


hV/UilO 


•IJtifcDE 
hi MITE 


•UtEDE 

hmiTE 


GEUGM 
AUPI No — ( Con 1 1 rim hJ 

F^nur 


h u IDKXi 


h PIUJIG 


hG'ivmj 


hilNDXC 


I- piJi'rtd 


r r.iF24 


hviEuix/r 


hVAlUO 

I 

— I*’IU.1N 
-PCAL 
-ilUMlD 

-lit-UIUM 

-G14-:UUM 
- t’Eb’llEE 
h t-YAUlO 
h lYAlUC 

h*VlNU2 

r WINU2 


■GmY4 
h LJ1’1)X4 

«.)U iK'Xt pago) 

— piYAKlO 


PGIG»«/\ 
h PGIGMA 
FIU.IN 


-L)KEl)E 

-UKITE 


■UtlDE 
h IXU'l'E 


ri]uCAT!<: 

— PsA'IXIIiS’ 

— pHvn'l{.N 


•UUU)E 

'MITE 


^KILLIN' 

F KILLIW 


H IVftil i'E 

-CUMPO 

-CUAlPO 


t . t . I 


-^A'FUHN 

-^UJCATI-; 







-KINUliX 

-ATtsru 

-WtliSlG 
-UPK2UU 
-KINUEX 
-ATliJ'fU 
-PKliblG 
-UPK2BU 
-KINUEX 
- ATliJ'IU 
-PUliSIG 
-ItliLUE# 

-SATSND 

-KINUEX 

-ATESrr 

-TEltP 


hWINUPrt 


AUPINS 


- INUPAlt 
- INUPAH 

- INUPAH 

- INUPAH 


oATUHN 

VHEUE 


-0ja 


CUGRIU 


-UIFFRS 

-UIFFHS 


LWADP 
hOIKADP 

•I'VAHIO 


NMCHIT 

Lhjyi 

-LEQl 

“DlFhHG 


.;''s 
• aa-Ai ' ' 

■ ; i 

-JGGAN 


-.4APP 

-GE'fl-TGW 


-NMCSOl 

-lQUS 

-IjQEU 

-JMOOTH 


-NMCHIT 

-QE 

- bVAHIO 


-UHIDE 

-DHITE 


rixT 

pGHIUCU 
[-GCUIGT 

pUUC ^CAIJi 


-NEX^r 
-GHIUCU 
-GCDIGT 
■GCUIH 


-SMUOC pSHUM 

— tJLJAnc’i 


SHAPFL 
f-SHAPFL 
r G3M009 

hUilEDE 

□KITE 
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b.a aubroatxno Uascription 


Tills section contains a dutatled description of each subroutine in ttie 
objective analysis program* Unless otherwise indicated the information is 
appropriate for tlie objective analysis witli the 4® x 6“ 4th order inodoi. 


SUBROUTINE ADPINS 


Argument list : Ncme 


Description: Perfortns obJt>ctlve analysis of sea level pressure and temperature 

and humidity, wind, and ge<^tential height at the 12 mandatory pressure 
levels; called by ALTER2* 

The routine also reads in all observatlcmal data, checks for time, report 
type, and quality marks, stores appropriate values in OQMMON, and calls respec- 
tive SCM analysis routines to modify model quantities. 

Restrictions: None 


Arguments and 
arrays passed 

Type 

Dimensions 

sT 

p2 

a3 

Description 

None 

Oormon areas 


Length 



Description 

CNTRL 



DBO 



4th order model parameters 

SATRCT 



54 



Satellite retrieval DSRN 

PHYSIC 



20 



Physical constants 

tJTYI'E 



4 



CJu^tity analyzed; 
set in ALTER2 

ADPTYP 



C 



Analysts parameters 

FAOOMM 


107800 



Analysis arrays 

ADPQWrr 


5DC4 



Data quality weights 

INDEXl 



320 



Report type descriptions 

INDF3(2 



64 



Report Type flags 

External references 




Description 


QE 

Indefinite logical conportson 

LEQl 

Single logical comparison 

TERP 

Horizontal bilinear interpolattori; 

PAIRZ 

Aircraft pressure calculation 

ATESTH 

Hydrostatic check 

ATESTP 

Sea level pressure check (bogus if necessary) 

ATESTt} 

Tests quality marks 

ATESTU 

Tests upper air Il-b quality marks 

ATESTT 

Tests TIR08-N quality marks 

ATEOTW 

Tests Wisconsin wind Il-b quality maiks 

ATEST3 

Category 3 report check 

FINDLV 

Variable pressure wind interpolation 

GBTAHT 

Aircraft height calculation 

INSEAL 

Sea level pressure SOI 

INSURF 

Sea level temperature SO! 

INUPAR 

URjer air SOI 

K INDEX 

Report index identifier 

PllESIG 

Mandatory pressure level sigma calculation 

RELDEW 

Relative humidity calculation 

UPKST3/UPKST5 

Il-a unpacker for surface/upper air 

UPK2BB/UPK2BU/UPK2BW 

Il-b unpacker for buoy/upper air/Wisconsin winds 
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SUBHCXITINh: AUPliMd (Continued) 


WlNOni Wind vector error check 

dATiM) Heads satellite sounding 


Iniait/Outtait ddiaune Pi I C 1 


20 X 

UJSAT(l) , . . . .LUSAT(NSAT) X 

81 X 

82 X 

83 X 

84 X 

8b X 

8(J X 


Method ~ ' liet^iption ~ 

U sei|. Il-a AUP data 
U soq. Satellite retrievals 
F seq. Drifting buoys 
F seq. Wisconsin Hast and West 
F seq. Wisconsin Indian Ocean 
F seq. Oonstant level balloons 
F seq. NAVAIUS 
F seq. Dropwitxlsondes 


‘ Stoivd 

2 

^ Fetcliod 
' Acidi’essod 


^UBKUU'l'INk: ALTh;K2 


Argument iiut. Mono 


f. 

I Uuscription: Driver routine for the 4tn orcScr objective analysis, called by 

^ MAIN, ftus routine interpolates tJie first quess to the pressutv levels, 

[. performs geostropluc adjustment if necessary, perfoitns objective analysis, 

[ interpolates analyzed fields to tne o levels, and ri^tores new values to tiie 

L model. ALTliKfi also performs dynamic initialization if necessary. 


Restrictions : None 


■} 

t Arguments and 

IVpu 

Diitiensions 

S 

F 

A 

Description 

arrays passed 








( 


Ooninon areas 

Length 

Description 

CN'ntL 

D50 

4tli order model parameters 

yfrYPli 

4 

Data flag set in the main loop 

AluPGr 

10 

Analysis control parameters 

CYCLb’ 

C 

Dynamic inialization parameters 

GDOSL 

4 

Geostropnic adjustment flag 

INS 

4 


INMAP 

4 

.vLip Index 

FGQW'r 

4 

First guess quality weiglits 




Description 


PAHA 

INUMN 

PUTWS 

AUPlNii 

GfcXJtSAD 

HEADIN 

rtINDXl 

KINDXii 

TWKITt: 

TAP24 


Initi.aiizes parameters and n,*ads ruuneiist 
Interpolates model o values to pressure levels 
Dynamic initialization restoring cycle 
Wsrforms objective analysis 
Geostropluc wind adjustment 
Reads mandatory level initial conditions 
Fetclies model quantities, reversing dbnetjsions 
Saves model quantities, reversing dimensions 
Writes model o values to tape 
24 tape 
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aiWHOITINb: ALTiiH2 (2.5* x 3* 2nd Order Model) 


Argument list: None 


u<^}cription. Driver routine for the 2.5* x 3* 2nd order objectlvo 
aiuUyeis, called by MAIN. 


Keistrlctions; None 


Arguments and 

Type 

Dimensions 

3 

F 

A 

Description 

arrays passed 








None 


Codinon areas 

Length 

Description 

dPLir 

24 

Split-grid flags 

AloPUT 

8 

Control flags 

GWOGl 

4 

Geos trophic adjustment flag 

TGAM3 

12441U0 

Model 0 quantities 

PHYGIC 

72 

Physical cotistants 

l/lYPE 

4 

Quantity flag 

duai’4d 

177GG 

Diagnostic quantities 


Llxt^enml references 


U t^cfiption 


tlDl, 

Uh'vV 

PAilA 

WtiTP 

WHTl 

INCvLN 

tiuauT 

raor 

ADPINa 

QjOCKa 

GhXJtSAU 

PTOGIG 

rtuiaK2 

IdiADIN 

'rwrtiTK 

UNSP 2 D 

UNSP 3 D 

WlJiaKi 


lieads model quantities from direct-acc«3ss disk 

Hewinds tape i 

Initiallings parameters for analysis 

Writes mandatory-level output to tape 

Writes model quantities to direct-access disk j 

Interpolates model o quantities to inandatory i 

levels i 

Heads diagnostic quantities from disk t 

Writes diagnostic quantities to disk 

Heads observational data and performs analysis 

CPU timer j 

Geos trophic adjusbnent of winds 

Interpolates analyzed fields back to a levels 

Heads teiporary disk 

Heads mandatory-level initial conditions 
Writes model o quantities to tape 
Unsplits model quantities 

Unsplits model quantities I 

Writes ternporary disk 


Input/dutput ddname 

1 I 1 0 Method 1 

Description 

2b 

X 

X 

U se»q. 

LbVHXi 2A data 

8 

X 

X 

U S(jq 

History tape 

33 

X 

X 

U seq. 

Teti^porary storage for diagnostic 
quantities 

21 

X 

X 

U seq. 

Mandatory- level ifilti-'.I conditions 

47 

X 

X 

U seq. 

Surface quiu!tj.tie*s 
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OP POOR QUALITV I 



SUBROUTINE ATESTA 


Argument list: NL. P, T, KB, * 


Description; Tests observational tenpomture profile for stability; sometimes 
called ADPINS. 

Restrictions: Profile runs fran surface to top of atmosphere. 


Arguments and 
armys passed 

Type 

Dimensions 

S 

1 

H 

Description 

M. 

1*4 



X 


Number of levels in profile 







(at least 2) 

P 

R*4 

NL 


X 


Observed pressures 

T 

R*4 

NL 


X 


Observed temperatures 

KB 

1*4 


X 

X 


Rejected profile counter 

* 






Rejected profile return 

Oontnon areas 



Length 


Descripition 


PHYSIC 20 Pljysical constants 









^jUdHOUTlNt:; ATi^ki'ni 


Argumunt U«t; NL, P, L, f, Kd, * 


Duecrlption. BsrfornB tiydrostatlc check on observational tentperature proi'ilo, 
culled by AOPlNd. 


liustrictions : 

Profile runs xrcro surface to top of atinospiiere. 

Arguments and 
arrays passed 

IVpe 

Uiinensions 

1 

1 

1 

Description 

NL 

1*4 



X 


Number of iovols in profile 
(at least 2) 

P 

H*4 

NL 


X 


Observed pressures 

1 

K*4 

NL 


X 


Observed heigiits 

T 

K*4 

NL 


X 


Observe ; teit4)eratures 

KB 

1*4 


X 

X 


dejected profiJ.e counter 

* 






dejected profile return 

Qxinibn areas 




Description 

PHYBIC 



:<2U 



Physical constants 
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iiUaaOUTlNh’ ATtiJTP 


Argutmmt Liiit: (JL), PSL 


Dusct'ipUon: Tosts sua level pi'es&iure u^iiist I'ln^t guobs tor pousvble 
unx>aeous eatrieti and be^^usses Uiein when necebsary, called ey AUPINd. 


Kesstrictiona; 

I'Jone 






Arguioenta and 

^ ■ 

Type 

Udnenaiuns 

d 

F 

A 

Description 

arraya paaaed 



i i 


i- - — . 


O) 

rt*4 

2 


X 

X 

blartli coordinates of report 

PdL 



X 

X 

X 

bea level pressure 

Ouinnon Areaa 


Lentsth 


Uescryition 

MOObl. 


7fc2i40 


Model first i;uess arrays 

CHTllL 



iJbU 


4th order model parameters 

liixieruoi Heforoncos 




— 

ii^ription 



Uanputes difference between obaorvation and model 


^-4 




iiUBttOUTINl:: ATiiki'IV) 


Artsutik^nt list. .^L, ML, X, Kti 


1 . 

\ 

I 


\ 


k 

I 

I 

! 


Detscription: T»it^ level ll-a data quality raarKs, called ADPINd. 

ties trict ions: k^lity marKs are returned uy UPtCSlb unpacaer. AUP data 

iilled witti UiK)UU. , if poor quality. 


Arguments and 

Tyne 

Dimensions 


F 

arrays passed 






Description 


ML 

1*4 



X 

Numuer of levels 

ML 

1*4 




Dimension of arrays 

Ttsyr 

K*4 

IC 


X 

Quality mark array 

X 

rt*4 

(.ML, 5) 

X 


Data array 

Kb 

1*4 

5 

X 

X 

Hejected report counter 

Qxnmon areas 



Length 

~Ui^criptVon 


None 


l^te rriaT references , 7 ' J ” 7. 1 ' UescfiptYon , 7_ I ” 

Lii)l iiingle logical canparison 


f 



itihfiffiiitirinfTiii 
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aUdl^OUTINt; ATiidTr 


Ar^unkitit Ust. * 


Uusicrlption. Tustts TIHLX>-N quality marks, filled oitlier operationally or 
interactively on McIDAii, called by AUPIN;:!. 

Host r ic tions . None 


/Vr^uinents and I iVpe Dimensions d F I A 
arra>s pissed | | j 


Description 


kTiiN 


Quality mark I’etia’n for pool 
quality 


SUdlUXlTlNk: ATLkj'fU 


Arguinunt list. NL, ML, KThiStf, X, KM 


Uuscrlptlon; Tests upper air LA.'vel 11-b quality uiarKs, called by AUFli>(M. 

Hestrictions: (Quality marKa aio as returned by UPK2BU unpacaer. 

Data tilled wltu U1I999., if poor quality. 


Arguments and 
arrays passed 

Type 

Dimensions 


F 

•• • 

A 

Description 

ML 

1*4 



X 


Number of levels 

ML 

1*4 





Dimension oi arrays 

KTfclST 

1*4 

(ML, 5) 


X 


(iiuailty mark array 

X 

H*4 

(ML,b) 

X 



Data array 

KB 

rt*4 


X 

X 


Hejected report counter 
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aUriauUTINH ATtaT# 


Artjurotjnt iitit; KWAL. Kd, * 


Ouscription; Tosts Wisconsin wind quality marks, caliud by AUPINS. 
tics trict ions, (^iity marks arc returned by UPK2BB unpackur. 


Arguments atid 

TVpe 

Dimensions 

3 

F 

A 

arrays pas.sed 







Description 


KWAL 1*4 X Quality mark 

KB 1*4 X X Hejected report counter 

* Hejected report return 
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iiUUKOUriNfc; ATiiil’J 


Argumeiit list; TcJUT, KAT, *' 


UusurLption. 'iVists C>itekSory A Iduntli'iur mrKs cor category os well us 
quality, called by AdPINl:». 

Kest net LOUS, (ijuaiity marks are ivturned by tue UPKdIb uiipacker. 


Arguments and 

IVlw 

Dbneiislons 

S 

F 

A 

arrays passed 







Description 


TlibT L*l 4 X Went 1 Her mai'k 

I’^T 1*4 X Qitugory returned =• 1. CATl 

=« 2 . CA'l*:i 

* Rejected report retui'n 

Ujiiiipn aivi^ lengt li ' ^ 

None 

i^ternui relei^erices iii^crljptl^j^ii 

Li‘XjL Dingle logical ojnparxsuii 
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UUtiiiUUTlNk: CUiHlU 


Argudtont list: CD, X 


Duuct'Lption: Computeti uartli coorditiat :s of tho 4th ordor isrld points, 

itestrictions; None 


Argiii.tt;)tits and 
arrays passed 

1 

'rype 

Uinvnsions 

3 

K 

A 

Description 

U) 

rt*4 

2 


X 


liarth coordinates 

X 

rt*4 

2 

X 



Grid points roturuud 
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tJUtitlUUTINt:: aiKAOF 


Arguiikint List. fiJ, KQ, L 


Uusurlptioa: Performs lateral "buddy" clwcK on rejected reports, callt?d by 

INSDAL. INSUKF, INUPAH, and tfltlUPU. 

Tills routine averages ail reports wlttun a 5* radius of the rejected 
reports. If there are at least 3 of these reports and the dlffeixsiice between 
the rejected report and the average differences is within the cutoU criteria, 
tlie report is accepted. 

Hestrictions: None 


Arguments and 

Type 

Uitnensions 

3 

F 

A 

Description 


arrays passed 









NS 

Kg 

L. 

1*4 

1*4 

1*4 


X 

X 

X 

Total number of stations 
(^lality flag 
l^essuiv level 

OjRinon areas 



Long til 

Description 

FAODMM 



10780) 

4th order analysis arrays 

ADPSC 



8 

Surface error criteria 

ADPUC 



bDO) 

Upper air error criteria 

OdSHNT 



120) 

(Served liMS storage area 

INDfcXi 



320 

Heport type description 

lixtvjrnal references 


Description 

LigL 



Single logical comparison 
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tiUBHOUTINli CUXXB 


Argument Ll-st; I 


Description. This routine allows tlie user to determine tlie time elai^sed in 
execution. 


Hestrictions. None 


Arguments and 

'Pype 

Dlrnensions 


F 

A 

arrays passed 







I IM 


Description 


'Phe amount of elapsed execution 
time. 
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iSUfttOUTlNE (JUTADP 


Arguinent list IOC,a),VALA.OIP,CUT,HMiS, It<M£i,l{MSX, ItiAlSX, 
FLA 

Duticrlption. Finds ttie difference between tne observation and tUe first 
guess and canpares the difference witti the cutoff criteria, called by If'fcJiiAL, 
INdUHF, INUPAK, and WlNUPrt. 

Uestrictiofis: A report flag is filled with .FALSE, if it is an allowable 

difference. It is filled with a "C" if it is not an allowable difference 
but within d times the allowable difference and with an "H" otnerwise. 






' 


Arguments and 

Type 

Dimensions 

3 

F 

A 

arrays passed 






K-X 

1*4 

2 


X 

U) 

l{*4 



X 

VALA 

rt*4 



X 

DIF 

K*4 



X 

cur 

H*4 



X 

UMd 

U*8 


X 

X 

IKMd 

1*4 


X 

X 

llMdX 

l{*8 


X 

X 

ruMd-x 

1*4 


X 

X 

FI A 

l.*L 


X 


Cannoji^ areas 



Length 


INDcXL SA) 


itescription 


Report type 
iiirth coordinates 
Reported data 
Uifi'eretice from motlel 
Error criteria 
HMS of accepted data 
Number of accepted reports 
HMS of rejected data 
Number of rejected rept^rts 
i{ep«jrl flag 


Description 

Report-type descriptions 


SUattOUTlNk: UE1>'1LE 


Arisumint list: lLXSt<N, MAXtikL. MAXL/m, mFb;. lAVKIi. 
ICK 


Description: Tliis routine pert’onns DEk'INE FILE functions, whicti must be done 
before using direct access 1/0 on a ddd dataset. 

itestrictions: None 


Arguments and Type Dimensions 
arrays passed 


Description 


IDSKN 

1*4 

Dataset reference number 

MAXHEC 

1*4 

Maxiinum number of records 
in tbe dataset 

MAXL'IH 

1*4 

MaxinHim length of a record 

ITYPE 

1*4 

'lype of variabiu 

lAVElt 

1*4 

Tile associated variable 

lOK 

1*4 

Variable returned 



aUHHOUTINt: UH>X4, DFDV4 


Argument list: J, JM, I. IM, F, D 


Uetjcription; Diroctlonul Uurivative of a two-iUinenslonal fioiU oy 4tn order 
dlfforunuiug, called oy Gi:)Qd«l. 

Hestr lotions. None 


Arguniuntij and 
arrays passed 


'IVpo 


Dimensions 



Description 


DruX4 

rt*4 


X 


X direction derivative returned 

DlDYd 

11*4 


X 


Y direction derivative ivturtieU 

J 

1*4 



X 

Y ^rid point 

JM 

1*4 



X 

Y grid maximum 

I 

1*4 



X 

X grid point 

IM 

1*4 



X 

X grid maximum 

F 

H*4 



X 

Horizontal fiolci 

D 

K*4 



X 

Di.stance between grid points 


56 


aUtittiJUTlNb: UIFbUS 


ArgURkMlt ilut. JM, IM, CD, D, F, DIF 


Uuticriptioi). Finds diftuivncu butwuen observation at arbitrary oartii 
coordiiiatos and a two^dununsional modei tieid using bilinear inlerpoLation. 

Kestr lotions: None 


Arguments and 
arrays passed 


Type 

DimensicitS 

S 

F 

A 


Description 

- 



„J 






JM 

IM 

O) 

D 

F 

DIF 


1*4 X 

1*4 X 

K*4 2 X 

K*4 X 

li*4 (JM.IM) X 

rt*4 X 


Latitudinal dimension 
Longitudinal dirnension 
blartii coordiiuitos of obsorvation 
Observational '/alue 
Horiaontal model field 
Difference returned 


v^iindti areas ^ . 11 . _ I Lengttt De scri ption 

None 


internal ref etvnces Deisc fipt'. • -n 

TfcHP Bilinear noriisontal interpolation ot model field 
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iiUlJtiiOUTIWl:; iMtDti 


At\;umunt ixst. Ul, liStii, JBCIK, IHbT 


U'iscrlption: D/A reud routine 

ties t net KM w. Mono 


ArguinentH and 
arrays passed 

Type 

Dimensions 

i 1 

O 

- - . 
b' 

A 

Description 

u) 

1*4 





Lot;icul unit number ot tiie 

liJUi 

1*4 





dei^iient numlAU’ 

JdUF 

1*4 





duller to puss data 

iKt;r 

1*4 





rteturn code 

vtSliJ 

1*4 







^UBiiUrriNE UKli’K 


Ansutnunt Uut: U), IdfcO. JdUF, Miiili. IKbT 


LXscrlption: D/A write routine 

Kestrictions : None 


Argiiments and 

IVpe 

Dimunsions 


F 

A 


arrays pasted 








Description 


UJ 

IM 

Logical unit number ol 

IdhO 

1*4 

degment nuntfjer 

JBUF 

1*4 

tkil’fer to puss data 

IKbT 

1*4 

Ueturn cudi; 

MdbO 

1*4 




SUiit^TINb UUMMyZ.litiAHT.rrtiUUA Artfumeiit Licit: 
OuicrLption. Duuiqy returns to ropiace provloutily 


None 

used routines. 


Kestr Lct Lons : None 


SUBROUTINE 


PILLIN 


Argument list: 


P 


Description: Computes o-edge pressures given the surface and tap pressure; 

called by SIGTOP. 

Restrictions: Profile is bottom to top. P(l) must be pre-loaded with sur- 

face pressure and P(10) must be pre-loaded with the tc^ pressure* 


Arguments and 
arrays [»ssed 

Type 




H 

Description 

P 


10 

X 

X 


o-edge pressures! 
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SUBROUTINE 


FIND 


Argument list: PA, P, X, KEV 


Description; Interpolates linear in log p a o- level profile to any desired 
pressure; called by SIGTOP. 

Restrictions: Set KEY«0 to avoid extrapolation. 


Arguments and 

Type 

Dimensions 

S 

F 

A 

arrays passed 







Description 


FIND 

R*4 


X 

Interpolated value returned 

PA 

R*4 


X 

Desired pressure 

P 

R*4 

9 

X 

o-level pressure 

X 

R*4 

9 

X 

o-level data 

KEY 

1*4 


X 

Extrapolation flag 


C^iiinon are as Length Deserrption ~ ~ 

None 

l^terrna l ref ererici^ Di^cripticm 

!/CATR Finds the o-level pressures which straddle a 

given p level 


SUBROUTINE 


PINDHM 


Argunent list: PA, P, H, KEV 


Description: Interpolates linear in log p a humidity o-ievel profile to any 

desired pressure; called by SIGTOP. 

Restrictions: Set KEY-1 to restrict relative humidity to between 0.01 and 1.0. 

Set KEYK) for specific humidity interpolation. 


Arguments and 
arrays passed 

Type 

^ 

Dimensions 

S 

P 

A 

Description 

PINDHM 

R*4 


X 



Interpolated humidity 

PA 

R*4 



X 


Deslr^ pressure* 

P 

R*4 

9 


X 


o-level pressures 

H 

R*4 

9 


X 


o-level humidity 

KEY 

R*4 



X 


Humidity unit's flag 


SUl3H(XJTrN^; 


FINDLV 


Argunent list: CD, NLEX', XI, X2, X3, X4, X5, KEY 


nescr .ption: Interpolates variable pressure Qategot^ (pilot balloon) wind 

reports to mndatory pressure lex^ls; called by ADPINS. If height is given 
rather than pressure, the model first guess height is used for Interpolation. 

Restrictions: Profile nust be ordered from bottom to top. Set KBY=1 if 

height Is given or set KEY»2 if pressure is given. 


Ai'guments and 
arrays ()assed 

Type 

Dimensions 

S 

F 

A 

Description 

CD 

R*4 

2 


X 

X 

Earth cooixilnates 

NLRV’ 

IM 



X 

X 

Number of levels in profile 

XI 

R*4 

20 

X 


X 

U wind returned 

x:> 

R*4 

20 


X 

X 

Height or pressure 

X3 

R*4 

20 

X 


X 

V wind returned 

\A 

R*4 

20 


X 

X 

U wind profile 

X5 

R*4 

20 


X 

X 

V wind profile 

KIA' 

IM 



X 


Height /pressure flog 

Oonmon areas 



Length 


tescription 




7E240 


Model first guess arrays 


rpy erences ~ ] Description 


TOiP IX^rfomK bilinear horl/x»ntal interpolation of first 

guess 

GR'l’WND FterfontB vertical interpolation of wind profile 






iiUtitiCXJTlNt:: 


WAKIO 


Argument liet; UJ, 10, LbYKL, JBUF 


Oescrlption. IViis routine perforoKi the direct access I/O u£ two-dimensional 
4" X 5* records. 

Ktistrictions: Tlie dataset must he KUOl'M^K, LK1 Uj* 372, UjUA>i>A aitd must be 
initialized by CALL Olil’ILE (LU, MSliJ, 1080 , 2, LP. IHbT), where 
MSHKJO* (number oi levels). Ilie number ol levels cannot exceed (35/30)* 
(tracKs in dataset). 


Arguments and 
arrays passed 

IVpe 

Dimensions 

3 

K 

• 

A 

Description 

LO 

1*4 




X 

Logical unit number 

lo 

1*4 



X 


I/O flag,"0 for read. 







=1 for write 

LLVbL 

1*4 



X 


4“ X 5“ recoixl number 

JBUK 

1*4 

1 

X 

X 

X 

4" X 5“ array 

Ccninon ai'eas 



Lengtli 


Description 

None 







External rel’erenctjs 

... . . ^ . 

Description 

UHLDE 


333 direct access ivad 

UJUTh; 


333 direct access write 

Input/Uutput ddname 

1 I 1 0 1 Method 1 


Liescription 

LU 


X X 333 DA 4° x 

5" record direct-access 


iJUBltOUTlNli GO>Il{ 


Ai*gumt}MV list; P, Q 


Uoscription. Tills routine detennlnes ttie direction from one point on ttie 
Bai'th iiiouti a great circle to unotner point on tno Earth. 'Oie direction 
is returned as an angle between 0® and 360® (0 represents north and 90 csist. 

Ifestr let ions; 'Fhe direction is returned as 0 ii the great circle -jatmot be 
determined. 


Argiunetjte and 

Type 

Oiinensiuns 


F 



A 

iXiscription 

arrays passed 








iJOJlll 

U*i 


X 

Direction returned 

P 

tl*4 

2 

X 

Origin point in eartli ccxirdinatos 


dM 

2 

X 

End point in earth ctxirdi nates 


- -- 

• 


l^crjjjtioji 

PHYSIC 



20 

Pliysicai Ciinstants 


auuKOUTiNt: ckjhs'I' 


Argutnent list. P, Q 


Description: Determines the distance from one point on the Lkrtii along a 

great circle to anottior point on the Earth. 

Restrictions. None 


Arguments and 

'rype 

Dimensions 


F 

A 

Description 

arrays passed 








GOHOT 

R*4 


X 

Distance returned 

P 

R*4 

2 

X 

Origin point in earth coordinates 


lt*4 

2 

X 

End point in earth coordinates 




1 

SU««3UTINE Gb)08Al) Argument list, ^k5ne 


Uescrlptlon. (Virfonrs a goostropnic adjustment on the first gues^ wind 
field based on an analysis of the geopotential height field, called by 
ALTtlK2. 


rtestrtctioas; T>ie first guess geopotential height field must be on direct 
access disK. 



Oomnon areas 

Lengtii 

Description 

MOULL 

7li240 

First guess winds and analy-ied iieights 

FAUOMM 

330U 

Work array for heiglit corrections 

External references 


Description 

tVAUIO 

i)ir<.jct -access 

1/0 ut.Uty 

GfXJdM 

Ckiostrophic adjustment at the mandatory levels 

^ 1. 

I 1 0 1 Method 1 IXjscription 

ZA 

X SSS DA 

Analysis dinxit-access file 


I 


i 
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origi>;al page ijs 

OP POOR QUALITY 


suBRarriNB 


GBOSM 


Argument list L, DZ, U, V 


Description: GEIOSM performs a geostrophlc adjustment of the first guess 

wind field. 

Restrictions: None 


Arguments and 

Type 

Dimensions 

S 

F 

■■ ■ 

A 

Description 

arrays passed 







L 

1*4 



X 


Mandatory level 

DZ 

R*4 

(46,72) 


X 

X 

First guess height differences 

U 

R*4 

(46.72) 

X 

X 


U-Component 

V 

R*4 

(46,72) 

X 

X 


V-OotnfX5nent 

Coninon areas 


Length 




Description 

PRINTl 


54 

Sample grldpolnts for printout 

GIDSl 


318 

Constants and constraints for geostrophlc 




adjustment 


PHYSIC 


4 

Physical constant 

External references 


D^crlption 


DFDX4 4th order X-dlrectlon derivative 

DFDY4 4th order y-dlrectlon derivative 
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i'niiafiiiiifti 


suBRarriNF. c.HTrcrr 


Argument list ; Ncme 


Description: Interpolates o-level AT fixan mid -mandatory level AT and adds o- 
level AT to first guess o-level temperature to obtain updated ten^jerature. 
Cbmputes potential temperature to perform a convective adjustanent and finally 
converts to temperature; cialled by PT081G. 

Res tr ict tons : None 



OoRtnon areas 

Length 

Description 

CNTT^L 

82350 

Model parameters 

MODRl. 

7R240 

Model arrays 

PHYSIC 

14 

Physical constants 

I’M ID 

2C 

Mid-mandatory level pressure 

PRMAND 

30 

F>ressure level constants 

Btternal references 

Description 


FHi.IN 

Computes o-edge pressures 


PSIGMA 

Deteimines mid-layer o pr^sure 
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SUBROUTINE 


OCTAHT 


Argument list: CD, ALT, PA, TZ, LM, ZZ 


Description: Creates pseudo height profile from an aircraft tempera tvjre 

report and the model first guess height field; called by ADPINS. The hei^t 
profile extends above the aircraft report only and is not computed below the 
aircraft level. 

Restrictions: The profile runs from top to bottom 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

F 

A 

D^cription 

CD 

R*4 

2 


X 

X 

Aircraft earth coordinates 

ALT 

R*4 





Altitude (not used) 

PA 

R*4 



X 

X 

Aircraft pressure 

TZ 

R*4 



X 


Aircraft temperature 

LM 

1*4 



X 


Lowest level desired 

ZZ 

R*4 

12 

X 



Height profile returned 

Carmon areas 


Length 



Description 

CNTRL 


D50 


Model parameters 

MODEL. 


7E240 


Model first guess arrays 

PRMAND 


30 


Mandatory level constants 

PRMANM 


2C 


Mid-mandatory level constants 

PHYSIC 


20 


PLiyslcal coastants 

External references 




Description 


TERP 


Bilinear horizontal interpolation of first guess 
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SUBROUTINE GETTajP 


Argument list: 7<P, TP 


Description: Oonpjtee mandatory level temperatures ; called by TAP24. 

Restrictions; None 


Argianents and 
arrays passed 

Type 

Dimensions 

m ^ 

s 

F 

A 

Description 

ZP 

R*4 

(46,72,12) 




Height field 

TP 

R*4 

(46,72,12) 




Tarperature field 

Oomnon areas 


Length 


D^ription 

PHYSIC 



20 



Wiysical constants 
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iiLriitOUTlNc; ub'l’KiW 


/VTt5Uini}Ut list. tCy, L, KJW 


lX‘scriptiotj. ObrrRiW roada In first ^uoss woiijnts f<jr a given quantity and 
Ujvei, called by INdJiAL, INdUUF, and IrJUPAK. 

Kustr let ions. None 


Arguments arki 
arrays passed 

'fyp-* 

r * 

Dimensions 

r 

o 

K 

A 

Description 

Kg 

I *4 



X 


Data type (i-U) 

L 

1*4 



X 


Level 

hUW 

H*4 

(46,7ii) 

X 


X 

Wo i gilts array 

Ounron area^i 





L^crip t iqn 


None 


tUternai qr«n^^^ Uescript ion 

bVAKIO Direct 'iccess I/U oi 4° x 5“ recoixls 

.L P. J. - 

(ii X SSd DA Dataset containing first guess weigiits 



SUBROUTINE GBIWID 


ArgMOlsnt list: imi, PA, PL, ML, XI, X3, X4, X5 


Deecrlptlgn: PerfomiB a llnoar vertical interpolation of a variable pressure 
wind [nrofile to nandat(»y levels; call«l PINUiV. 


{Restrictions: None 


Arguments and 
^rays passed 

Type 

Dimensicme 

S 

F 

A 

Description 

NLEV 

1*4 



X 


lAjmber of input levels 

PA 

RM 

20 


X 


Input pressure 

PL 

R*4 

12 


X 


Output pressure 

NL 

1*4 

20 




Not used 

XI 

R*4 

20 

X 



Output u wind 

X3 

R*4 

20 

X 



Output V wind 

X4 

R*4 

20 


X 


Input u wind 

X5 

R*4 

20 


X 


Input V wind 
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oaiJ2-lAL PAGE Ifc 
OF POOR quality 


SUBROUTINE GRIDCD 


Argument list; JI, X 


Description: Converts a 4th order grid point to earth coordinates 


Restrictions; None 


Arguments and 
arrays passed 

Type 

Dimensions 

- - ........ 

' 

S 

F 

A 

Description 

JI 

1*4 

2 


X 


Grid point 

X 

R*4 

2 

X 



Earth coordinates 


SUBROUTINE 


GTOPOG 


Argument list : IDSRN 


Description: Reads topography values which are different from those used 

by the model; called by ALTOG. 

Restrictions: Vfhen the variable IDSRN in the argument list is less than or 

equal to zero, a return is eicecuted and the read is bypassed. 


Arguments and 
arrays passed 

" — 

Type 

Dimensions 

. i 

S 

F 

r* — — 

A 

Description 

IDSRN 

. 1*4 



F 


Data Set Reference Number 

OoRinm aiceas 


Lengtii 


Description 

7.KDEL 


7E240 


Model first guess arrays 

External references 





Description 


None 


Input/Output ddiuune 

"T" 

0 ! Method ' 

Description 

IDSRN 

X 

Seq. 

Sequential 
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page it' 
OP Pr>'H QUALITV’ 


SimUaiTlNK GTOF 


Argument list: 


None 


l>escrLptlon: Computes mid-mandatory level AT; called by PTO8IG. 

Restrictions: None 


Arguments and 

Type 

Dimensions 

S 

F 

A 

Description 

arrays passcxi 








None 


Corroon areas 


Length 


Description 


CNTKL 

MODEL 

PHYSIC 

PRMAND 


DBO Model parameters 

7E240 Model analyzed fields 

20 Ptiysical constants 

60 Mandatory level constants 
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SUBROUTINE HJTEjfi* 


Argument list: TS, P, Z 


Dtecriptlon: GonoxiteS mftnd&tot*jr level geopbtential hei^t given the 1000 tito 
height, mandatory-level pressure and mid-mandatory level temperature. 

Restrictions: None 


Arguiaents and 
arrays passed 

* 

Type 

Dimensicms 

S 

p 


Description 

TS 

R*4 

11 


X 


Mid-tnandatoiy level tefrqjeratures 

P 

R*4 

12 


X 


Mandatory level pressure 

Z 

R*4 

12 

X 



Mandatory level height 

Oeenibn Areas 


Length 



Descriptidi 

PHYSIC 


20 



Physical constants 
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SUBROUTINE HUM IP 


Argument list: PA, P, H, KEY 


Description: Performs a linear in log p interpolation of a mandatoiy-level 

relative humidity profile to the o-leveis; called by PTOSIG. 

Restrictions: Set KEY=1 to restrict the relative humidity between 0.01 and 1.0. 
Set KEY=2 for interpolation of ARH. 


Arguments and 
arrays passed 

TVpe 

Dimensions 

S 

F 

A 

Description 

HUMID 

R*4 


X 



Interpolated humidity 

PA 

R*4 

12 


X 

X 

Fissure profile 

P 

R*4 



X 

X 

o-level pressure 

H 

R*4 

12 


X 


Humidity profile 

KEY 

1*4 



X 


Humidity flag 


SUBROUTINE HON 


Argument list : Nt^e 


Description: Fills model arrays for analysis by vertical interpolation 

(using SIGTCX>) from 4th order o-levels; called ALTER2. 

Restrictions: Array T (J,I»12) is filled with 1000 mb height. 


Arguments and 

TVpe 

Dimensl(X)s 

S 

P 

A 

— 

Description 

arrays passed 








None 


OOinnon areas 

Length 

Description 

CNTRL 

D50 

Model parameters 

MODEL 

7E240 

Model arrays 

QTYPE 

4 

Data flag from ALTER2 

CYCLE 

C 

Restoring cycle parameters 

External references 


Description 

SIGTOP 

Single profile vertical interpolation 


SUBlttXJTINK INSEAL 


Arguniont list; NS 


Description: Sea level pressure SCM analysis; called by ADPINS. 

Restrictions : None 


Arg>.«nents and 
arrays passed 

TVpe 

Dimensions 

S 

F 

1 

Description 

NS 

IM 



X 


Number of stations 


(jommon areas 

Length 

Description 

PHYSIC 

28 

F’I'iysical constants 

ADl^TYP 

18 

Analysis parameters 

SC3VI2 

18 

SCM parameters 

CNTRI. 

D50 

Model parameters 

PRINTl 

54 


MODEL 

9B40 

Model arrays 

FACXWM 

140300 

Analysis arrays 

ADPSC 

8 

Sea level error limits 

SMOOTS 

18 

Shoothing parameters 

OBSKNT 

1200 

RMS storage area 

FTfQWT 

33C4 

First guess weights 

ADFQWT 

5DC4 

ADP quality weights 

INDRXl 

320 

Data ty|)e 


External references 


Description 


LBQl 

IQW 

LQUS 

SvSCAN 

aiKADP 

niTADP 

niFFRS 

FVARIO 

Nr«csai 

SMOOTH 


Single logical comparison 
hirope grid point determinant 
United States grid point determinant 
Station counting scanner 
Horizontal buddy check 

Check for observation and moool difference 
Difference between observation and model 
Direct access I/O for 4° x 5° records 
SCM accumulating scanner 
4° X 5° field snoother 


Input/CXitput ddname 


I 


0 


Method 7 


Description 


23 


X SSS DA Temporary direct-access dataset 
for the first guess and analysis 























SimROUTINE INUPAR 


Argument list: NS, KQ 


Df?scription: Upper air SCM analysis; called by AlPINS. 

Restrictions: None 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

F 

A 

Description 

NS 

1*4 



X 


Number of stations 

KO 

1*4 



X 


C^antity type 

Garrnon amas 


Length 


Description 

(S.'ime as INSRAL) 

External references 





Description 


QE 

LHJl 

rxjEU 

LQUvS 

SSCAN 

ClfiCAIP 

CITTADP 

niFl'liS 

f'VARIO 

NMCHIT 

meson 

SMOOTH 


Indefinite logical comparison 
Single logical comparison 
Europe grid point determinant 
United States grid point determinant 
Station counting scanner 
Hori/x)ntal buddy check 

Qiock for oJxservatlon and nvxlol dlffennice 
Ccxnputes difference l>etwc>en ol>>ervation and 
first guess 

Direct acc(«s I/O for 4° x 5° records 
Geopotential height stability cli€?ck 
SCM accumulating scanner 
4° X 5° field frnother 


Input/Output ddname 


I 


0 


Method 


Descriptio n 




X SSS DA Temporary direct-access dataset 
for the first guess and analysis 


S3 




i^UUHOITINt; KINDclX, UINDiiX Al^Ufflytlt Ust: CAT. fcJTATN, IWI. LW, OJ 


Ufi-scription. Heport type identifier, caXied by AUPINb’, 

flue routine tags a report witn an index number detenninod tra.i tin) lata 
category, station IlJ, data source index, ifistrument type, and earth ccxin • 
nates. tiicU separate report type corresponds U> a dirierent set of qm 
Weights. 


ties tr let ions: A aixj is returned if the report is unrecognizabie. 


Argutnents and 

h — - — — - 

'iVpu 

Dimensions 

r - - - 

a 

r ‘ * 

F 

A 

Description 

arrays passed 







MNUfcX 

1*4 


X 


Report index 

CAT 

L*i 



X 

category 

oTATO 

L*i 

6 


X 

atation ID 

ijai 

1*4 



X 

Data source index 

iNa 

1*4 



X 

Instrument type 

CD 

H*4 

2 


X 

bkrtn cootxlinates 

UiNDfcX 

K*8 


X 


Report description 


— 


lengtii 

Lfesoriptiq^^ 

None 






kixternal references 






Lfc)yi aingie logical c«>nparison 








iJUbliOUVlNli 


Ari^uiUMit UvSt. A, tj 


l isjci'iption. 0><iiparots two aitigle byte cuaraotory, U' they arc tiic sajnc, 
iticti ftiturtia. TKUi'A, oUi-jtvi.se tt f-,-turas. 

Hes t r Lo 1 1 otus : Ncmic 


Arguments arui 
arrays passed 

IVpii 

1 

Du.'KJilsioris 

o 

F 

A 

1 Jose ript ion 

bli^i 

1*4 


X 



IjOglciil comparison 

A 

i.*L 



X 


ijliur;ti;t-,'r L 

li 

b*l 



X 


Character 2 


-So 


SUBROUTINE LOC 


Argument list; MAP, II, JJ, NX 


DeRcrlption: Selects the quantity to be plotted by testing the logical array 

QMAP; called by MAPP. 


Restrictions: None 


Arguments and 

Type 

Dimensions 

S 

F 

A 

Description 


arrays passed 









MAP 1*4 (46,72) 

II 1*4 40 

JJ 1*4 72 

NX 1*4 


CbRincm areas 


QMP 


External references 


SCALE 

PRSCAL 


length Description 

50 Maximum number of maps 


Description 

Plots the quantity on a map 
Prints >he scaled map 
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0 " 







SUBROUTINE l/)CATE 


Argument list: PA, P, N 


nescrlptlon: Finds the pressures which straddle a given pressure level. 

Ijov*;! returned is the level above the given pressure. 

Restrictions; Profile must be ordered from bottom to top. 


Arguments and 
arrays |xisse<.l 
1 

Type 

Dimensions 

3 

F 

— 

A 

Description 

'iCX:ATE 

1*4 


X 

X 


Level returned 

PA 

R*4 



X 


Given pressure 

P 

R*4 

M 


X 


Pressure profile 

N 

1*4 



X 


Number of levels 
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SUBROUTINE LQEU 


Argument list: J, I 


Description: Determines whether 4** x 5** grid point Is over Europe. Europe is 

defined as the region from 38®N to 66“N and ’ran 5*W to 30® E. IjQEU is .TRUE. 
If the grid point falls In this region. 

Restrictions : None 


Arguments and 
arrays passed 


Type I Dlmensl<xis 


D^crlption 


Logical flog returned 
latitudinal grid point 
Longitudinal grid point 









i 


SUBROUTINE LQUS 


Argument list: J, I 


Description: Determines If a 4* x 5* grid point is over the United States. 
The United States is defined as the region from 26*’N to 58®N and from 135“W 
to 80®W. liQUS Is .TOUE, if the grid point falls in this region. 

Restrictions: Ncmio 


Arguments and 

Type 

Dimensions 

S 

D 

1 

Description 

arrays passed 




■ 

■ 




I 

'i 


IQUS 

J 

I 


L*4 

1*4 

1*4 


X 


X 

X 


Logical fl<>g returned 
Latitudinal grid point 
Longitudinal grid point 






iiUBKiJUTINI:: MAPP 


ArifUfnent list; KNT 


Usscriptl<¥i; Increments an integer variable whiub dticiUes the pressure level 
at ’I'icb tbe map is to printed, 'this is called by INdEAL, INiSUlit'' and INUPAH. 


t.r actions; None 


Arguments and 

Type 

Dimensions 

d 

p 

A 

Description 

passed arrays 








KNT (46,72) Araray of integers 

Oannon ^ rea s rTT-II LengU i ix^ription 


U«iSAP 

I^ei^r refereitk^ 


ux; 


E.llg 


ri ettMXi r 


Alap number 

PescHiption 

delects tbe map to be plotted 

ET- - I - 


SUBROUTINE NESITTPE 


Argument list ; NT2 


Description: Sends message to curator to mount a new tape; called by TAP24 

Restrictions : None 


Arguments and 

Type 

Dimensions 

S 

F 

A 

Description 

arrays passed 








NT2 R*4 Logical unit nuiriber 
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SUBROUTINE NE9CT 


A3PgUTV?nt list; DMAX, CD, I, J, * 


Description: NEXT is the grid point scanning routine; called by SSCAN and 

NMCSC3/!* This routine is used to efficiently scan all grid points within a 
given radius about a point on the Earth. NE9?r is called repeatedly returning 
a single grid point each time until all points are read in the given radius. 

Restrictl(Xis: Beware that conventional GLAS definitions of I and J are re- 

versed. Also, all the points in a given radius should be read before moving 
to the next point. 


Arguments and 
arrays passed 

TVpe 

Dinvj'^sions 

S 

F 

A 

Description 

DMAX 

R*4 



X 


Scanning radius in km 

CD 

R*4 

2 


X 

X 

Origin in earth coordinates 

I 

IM 


X 

X 


Latitudinal grid point returned 

J 

1*4 


X 

X 


Longitudinal grid point returned 

lit 






Return when scanning is coi^lete 

(Icnitxxi areas 


bsngth 



Description 

BNDS 



8 



latitudinal hounds 

CJfTRL 


D50 



Model parameters 

PWSIC 


24 



Physical constants 

Elxternal References 





Inscription 

DOR ID 


(3onverts earth coordinates to grid points 


SUBROUTINE ^MCHIT 


Argument list: NLEV, NITER 


Description: Oomputos the model first guess hel^t field; called by INUPAR. 

Restrictions: Ehch level must be called in the order given in INUPAR. 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

F 

A 

Description 

NLEV 

1*4 



X 


Pressure level 

NITER 

1*4 



X 


SCM iteration number 


OcjtTmon areas 


liei^tlh 


Description 


CNTRL 

FAOOMM 

MODEL 


D60 

107800 

7E240 


Model parameters 
Analysis arrays 
Model arrays 


■ * 

! 

ittiil liili II It! I I iW ■ 
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SJUBHOUTINr) i^vlLkiOI 


Arsuinent List: JM. IM, Him, HADPAC, OJ, DIP. ^T. 
WDIH. Dhai, CN'r, HAD 


Description. Huccossive corrections accumulator, called uy INSEAL, INdUHP, 
and INUPAH. 'Alls routine accumulates tlie point corrections, apply itig 
tUe appropriate distance and quality weiglits. 'Aie wind correction is also 
dependent on direction with upitreain and downstream ousurvations given more 
welgtit titan crosswlnd observations. Aie scanning radius is decided on density 
ol stations about a grid point. 

Kestrlctlons: None 


Arguments and 
arrays passed 

Type 

Dimensions 

d 

i - - 

F 

A 

Description 

JM 

1*4 



X 


Latitude dimension 

IM 

1*4 



X 


Longitude dimension 

Him 

H*4 



X 

X 

Maximum iteration radius 

HADPAC 

H*4 



X 


Iteration distance factor 

CD 

H*4 



X 

X 

Hartli coordinates of observation 

DIF 

H*4 



X 


Observational and first guess 







difference 


H*4 



X 


Observational quality weight 

WDIH 

H*4 



X 


Obseivatlonal wind direction 

DHL 

H*4 

(JM,IM) 

X 

X 


Corrections array 

iwr 

H*4 

(JM.L4) 

X 

X 


Weights array 

HAD 

H*4 

(JM,IM) 


X 


Hadii array 

Cbnnion areas 

■ ~ ' • — 

Length' 



Da^r^i^tion 

PHYSIC 


20 



Wiysicai constauLs 

Ext^nal 'i^ferenc^ 




Description 


Grid point scanner 

GCDIH Great circle direction calculation for elliptical 

weigiiting 

GCUIiff Great circle distance calculat ion 

GKIDID Converts grid points to earth coordinates 




SUBROUTINE PAIRZ 


Argument list; AZ 


nescrlption; Computes aircraft pn^ure from reported pressure altitude 
and U.S. standard heights; called ADPINS. The height Interpolation 
Is linear In log p* 

Res trictlons : ?tone 


Arguments wid 

Type 

Dimensions 

S 

F 

A 

Description 

arrays passed 








PAIRZ 

AZ 


R*4 

R*4 


X 


X 


Pressure returned 
Reported pressui*e altitude 



SUBBOUTINE PARA 


Argianent list ; None 


D^criptlon: Initializes the analysis routines; called by ALTER2. 

Restrictions; None 


Arguments and 

Type 

Dimensions 


D 

A 

Description 

arrays passed 




■ 




None 


Ocranofi areas 



De«criptic«i 

CNTRL 


D60 

Model parameters 

Auvr 


C 

ATLER2 parameters 

SATRET 


2C 

Satellite retrieval datasets 

1NDEX2 


64 

Report-type flags 

PHYSIC 


24 

Physical constants 

ADPTYP 


10 

Analysis parameters 

StSClOTS 


18 

Surface smoothing flags 

SMOOIU 


2D0 

Upper-air smoothing flags 

PRMAND 


90 

Pressure level constants 

PRMANU 


90 

Reversed pressure level cOTstants 

PRMANM 


90 

Mid-pressure level constants 

ADP8C 


8 

Surface data cutoffs 

ADPUC 


5DC0 

Upper-air data cutoffs 

80(1 


. 18 

parameters 

GB081 


318 

Geostrophic adjustment parameters 

PAOOM 


170 

Work area 

ADPQ^ 


5D04 

FGGE data quality weights 

INDEXl 


320 

Report-type descriptions 

PRINTl 


54 

Grid points printed out 

External references 



Description 

DINDEX 


Report type description identifier 



Description 

5 

X 

NL seq. 

ALPUTZ namelist 

26 

X 

U seq. 

RMS output 

23 


X SSS DA 

Temporary direct access dataset 
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suBROurrMT-: k’Ai. 


Argument list: P, PBAR 


TV?scrlptlon: Computes mid cj-level pressures frcm o-edge pressures* 

Restrict ions : None 


Arguments tuid 

Type 



Dimensions 

S 

F 

A 

Description 

arrays passcxl 








X 


P 

!>BAR 


R*4 

R*4 


10 

9 


X 


o-level edge pressure 
Mid o-lovel pressure 


SUBKXfriNE POTBUP 


p 


1 


Argument list: P, TBAR, TS 


Descript Icm: Oonputes potential temperature at mldHnandatory levels • 


Restrictions: flone 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

1 

1 

Description 

P 

R*4 

12 


X 


Mandatory level pressure 

TBAR 

R*4 

11 


X 


Mid-mandatory temperature 

TS 

R*4 

11 

X 



Potential tam»rature 






SUBROUTINE PRESIG 


Argument list : P 


Description: Determines mandatory level number from pressure; called by 

ADPINS. 


Restrictions: Msindatory levels are numbered from top to bottom. 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

F 

A 

Description 

PRF23IG 

R*4 


X 



Mandatory level pressure 

P 

R*4 



X 


Pressure 


SUBHCUTINE PSIGMA 


Areument list: P2. PO, P 


Deecrlption: Computoe the mid-level pressure of the o layers fran the 

surface preseure and the top pressure. 

Restrictions: P(10) is returned as surface preamire. 


Argumsnts and 
arrays passed 

Type 

Dimensions 

S 

P 

A 

Hescrlptlon 

P2 

R*4 





Surface pressure 

PO 

R*4 





pressure 

P 

R*4 

10 




<3 pressure 


SUWKXITINK I’SURFE 


Argument list: Mono 


I\>scrlptlon: Computtjs surface pressure from sea level pressure, sea level 

tttnperatuix*, and surface height ossunlng a standard lapse rate; callal by 

Rc«?trtctions: Mono 









Ar\;uments and 
/"rays Passed 

Tj'pe 

Dimensions 

S 

F 

A 

Description 


Ncxh? 


Oomlon Areas 


Length 


Description 


CKTRL 

M(T)ra. 


reo 

0!MO 

20 


Model parameters 
Model arrays 
Physical constants 


SLTBRarriNR PTOPIG 


ArBvuni*nt list: None 


rk>sorlption: Intcrfx>latoR AT, ARH, Au, anrt av to tlie o levels and odds them to 

first gu««s fields; called by ALTFIR2. 


Restrictions: 

None 




An»umi’nts and 

1 

Type 

Dimensions 

S 

F 

armys passed 







Description 


None 


Cdnnon areas 

Length 

Description 

cjraii. 

82050 

Model parametere 

MODEL 

7E240 

Model arrays (P and o) 

QTYPE 

4 

Data flag from AIjTER2 

PRMANU 

30 

Mandatory level constants 

roowT 

33C4 

IXimny array 

Bctemal references 


Description 


PCAL 

Computes mid-level o pressure 

GTTOT 

Computes mid-mandatory AT 

HUMID 

Interpolates relative humidity to o levels 

WIND2 

Interpolates wind to o levels 

FIIJ.IN 

Computes o-cdge prossut\‘ 

Gixyror 

Interpolates AT to o levels 

FVARIO 

Reads 4° x 5® data records 

PSURFE 

Computes surface pressure from sea level pressure 

SPFiaiM 

Converts relative humidity to specific humidity 

GFxmir 

Interpolates AT to o levels 


dUHrtOUTiNfc; purwd 


ATtiuinunt 1 ist : 


LXiscrtption . l^tiaj.uc initlait^ation restoring cycio, cjailed by Al/riiH 2 . 

'ft;' routine maxes sucoissive caiis to tue incxiel uyarotlytimnic« witJi 
alteriiat „iy forward and bactoward tune ste^JS tn order to balance ttie masjj and 
wind r'Letds. 


K« -‘Strict ions. 

None 




Arguments and 

'iVjJe 

Ddniinsions 

8 

F 

- 

A 

r*"*‘*“‘“*' 

Description 

arrays ptissed 



^ 





IM 


X 

Ai;riiK2 data flag 

Cufif.ion areas 


Lengtti 

Description 


bUJHTll/dVin’iVyMOm/iWiVKIV.SiKi ii</W(jiK/ttfUiKii/ 4tb order model 

diviuo/ i’wj/.vu>y/CN'i’dJ. 


DllUO- 

CYCUi 


4 

b7540 


C 


Debug flag 

Nkxiel analysis array.s 
Hestoring cycle flags 


lixLurnal rcierencus 


(JOVIIJU 

TV^'UTl'J 




/todel liydrodynajiiics 
Uutijuf. U) tj(JLi8 tape 


103 


dUtitUJUTINI:: 


ArguiiKjnt list. /A, 62, NB 


Uuscrlptton. Campdrison ol' two cUaractor strlrtgs of a liivon lou^ith. If ttio 
two atrings at>* ulontui-il, roturna .TBUH., otiioiviso vit; I’otunia .KALiih’. 

itistrLctiono; Nono 


Ai'^umonts and 
arrays passoU 

Type 

Uiinunsions 



A 

Description 

Vjli 

L*4 


X 



Uogicai canparison 

ZL 

L*i 

NB 


X 


Oiaractor strin<i 1 

Z2 

b*i 

NB 


X 


Oiaracter strintj 2 

NB 

1*4 



X 


ixjngtli of cliarajtor stritigs 
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► 





i 


iiUHJiOUTLNK ttiJ-VDIW Artjuinetit list. NTJ, TAUI, N 

Uoscriptioii. Heads lu mandatory- level initial conditions for analysis, 
culled by AUTHH2. 

I lust net tons: i’Jie datusot ftiust ixj tn tne standard nuiKiatory-ievel format. 


IrgumutiUj and 
arrays parsed 



'iVpu 

Dua^nsions 

li 

K 

f 

A 

lX»scription 

N'V2 

IM 



X 


DSHN for initial c<tnditions 

TAU 1 

li*4 





(not used) 

d 

1*4 



X 


nata fluti from /UjTDHli. 

Common a^£W 

- - 

Lengtli 

- ■ 

- - - 

Description 

CNTHL 


UbO 



Model parajfleters 

MUDKlj 


71i2^lO 



Ivlodel array.s 


reljuivnc^^^^ , . . T- . ] . ' LX^criptxon 


I Noiii* 

^ (^t£Ut ddnujiK} , , ] I .[ O . . ' ” ' V. ] . 

N’I‘J 


X 


U si*q. 


Muiidulory-level itiittul coiidi tiotts 


SUBROUTINE RELDEV 


Argifflient Hat: P, T, DPD 


Description: CbnqHites relative humidity from pressure, tenq:erature, anu dew 
point depression; called by ADPINS. Relative humidity is returned as a value 
between 0.0 and 1.0. 


Restrictions: 

REIDE^ is filled with 09999. 

, if 

invalid arguments are pm^J. 

Arguments and 
arrays passed 

Type 

Diitienslhns 

i' j I’ 

i 

A 

Description 

’ „ r - - -,r ■ r ^ - . 

RI3DEW 

R*4 


X 


Relative humidity returned 

P 

R*4 


X 

X 

Pressure 

T 

R*4 


X 

X 

Tanperature 

DPD 

R*4 


X 


Dew point depression 

Oenroem ar^is 



length 


Description 

None 

External References 



Descript icxi 


SATURN SaturaticHi specific humidity I 


SUBROUTINE RELHUM 




Argument list: H, T, P, L 


a specific humidity profile to relative humidity 


Dimensions 


S 


Description 


L 

L 

L 

L 


X X 
X X 
X 


Iftjinidivy profile 
Temperature profile 
Pressure profile 
N>jntoer of levels 


Length 


Description 


Description 

Saturation specific humidity 
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SUBROUTINE RINDXl Argument list: NQ 


Description: Pills analysis arrays with nodel arrays, reversing the grid 

Indices. 


Restrictions: None 


Arguments and 
arrays passed 

Type 

^ 1 . „ „ 

Dimensions 

s 

1 

A 

Description 

NO 

1*4 



X 


Quantity flag ftxm AIjTERS. 

OcinncNi areajs 



HHi 


Description 


CNTRL R3!)30 ^fodel arrays 

MODEL 57550 Analysis arrays 










SUmUXJTINE liINa)C2 


Argument list. 


Description 

indices. 

Kcstrictlons. 


Arguments and 
arrays passed 


Odiinnn a reas 

CNniL 

W0KK2 

MaitL 

CYCLE 


Fills model arrays witli analysis arrays reversing the grid 


None 


^ 

— — . . • . « - 


r • 


^ «... * ^ 

Type 

Dimensions 

S 

P 

A 

Description 


1*4 X Quality flog fran Al..TEa2 

Length lii^cription~ 


83D30 Model arrays 

81600 Scaled model arrays 

57540 Analysis arrays 

C Kcstorlng parameters 


ion 



SUBiW:lNU 


ArgvinBnt liat KU, LU, TAU, m’lU, D8I. 

KINS. ZLAT, ZliQr^ KTSN, 
ai!V. P3K, TSF, DPCSF, 
1>SL. ZNTH Ur. MK P-fP, 
rrp, Tun, WLP, KW ZM, 
TM, DRUM, • 


Ooscrlption Reads satellite sounding reports called by ADPINS. 
Restrictions : None 




' 





Arguraeiits and 
arrays passed 

lype 

Dimensions 

S 

1 

F 

A 

Description 

KU 

1*4 



X 


Format ID * 1,2; or 3 

UJ 

1*4 



X 


Uogicai unit number 

TAU 

R*4 


X 



Hour of retrieval 

STNID 

R*4 

2 

X 



Statiot» ID = blank it missing 

USl 

1*4 


X 



Data Source Index 

KINS 

1*4 


X 



Instrument type 

ZIAT 

R*4 


X 



Retrieval latitude (-90" to 90") 

ZLON 

R*4 


X 



Retrieval longitude (0“ - 3U0“W) 

KTSN 

1*4 


X 



Quality mark 

tmv 

R*4 


X 



Station elevation 

PSP 

R*4 


X 



Surface pressure 

TSF 

R*4 


X 



Surface temperature 

DPDSF 

ft*4 


X 



Surface dew point depn?ssion 

PSL 

R*4 


X 



Sea level pressunj 

ZNTH 

R*4 


X 



Zenitli atigle 

IOC 

1*4 


X 



tliannel code 

MR 

1*4 


X 



Retrieval code 

p^rp 

R*4 


X 



iropopause pressure 

TTP 

8*4 


X 



Tropopause temperature 

'IlM 

R*4 

15 

X 



Mean layer temperature 

WLP 

R*4 

3 

X 



Precipitable water 

KTW 

1*4 

18 

X 



Temperature and precipitable 
water quality marks 

ZM 

R*4 

15 

X 



Mandatory level heiglits 

TM 

8*4 

16 

X 



Mandatory level temperature 

DPDM 

8*4 

6 

X 



Mandatory level dew point 
depressions 

Osnniou anreas 

" 

‘Leng’tb 


' ~ 1 1 

None 







External references 

- 


— 


Description 

VHIDE 





Variable reading utilty 



SUBROUTINE SATSND (Continued) 


Input/Output ddnawe 1 I [ 6 1 Method T 
III X U Seq. 


Description 


TIROS-N database 


\ 

I 

) 

I 


I 


f 


iv 

•; : 

i-f** 






Ill 




SUBROUTINE SAtVHN 


Aivunuiit Uat. T, P 


i Uuscrlptlon Canpatus saturation spocU’ic hunidity in i$/k Cixin tu>npurai:ur*u 

[. unci pn^uro. 


IU3Strictluns. Erronixais lvalues may causo prugrain chocks. 


Arguments and 
arrays passed 

H 

Diinetisioits 

r— ■■ 

S 

F 

A 

fXjscrlption 

SAIl/RN 

ft*4 


X 



Saturation specific huioiclity 

T 

H*4 



X 


'ftjniporatatv 

P 

R*4 



X 


Pressure 







sumttJUTiNt: iBLizii 


Argunuut list. TS. P 


IXascription Stabiliziis the potential tonvwraturu profile at the mandatoi*y 
levels. 

Kestr let ions; Profile stiouLd be oixlered from bottom to top. 

Ari^uments and 
arrays passed 

Type 

Dimensions 

L 

s 

F 

r-*” - 

A 

Description 

T3 

«*4 

11 

X 

X 


Itotential tetqperaturt; 

P 

lt*4 

12 


X 


Pressure profile 

Cannon areas 

" 

Length 

- - 


Description 

OOUOT 


S 

1 



Unstable profile counters 

External references 




Description 

coNsrr 



Stabilizes 

unstable profile 
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SUBROUTINE SCALE 


ArfUBirat list: MAP, IR, JR. NX 


Deecrlption: Plots th© al^uibetlc cteraeter or dicit depending on the 

frequency of th© quantity and the choice of th© maxlmuin number of levels; 
called by UC. 

Restrictions: None 


Argunent and 

Type 

Dimensions 

S 

P 

A 

Description 

arrays passed 








1*4 (40,72) 

1*4 Latitude Index 

1*4 liongltude Index 

1*4 Meteorological quantity to 

be printed 


114 







SUUlttXJTlNE SJIAPFL 


Ar^unuiit List F, NSM TWIND 


Lkiscriptiotj Shaplixj filter applied to a horikioiital field called oy S&txrni 
No :smootlimg le done un tUe pule points. 


Uestrictions. Set TWINI) to .TltUE. if aiuxitliintj tlie wind field. 



r - - - - - “ 




r " * “ r 

Ar^unents and 

Type 

DDnensions 

S 

F 

A 

arrays passed 







Description 


F 

NSM 

rwiNi) 


KM (46,72) X X 
IM X 

L*4 X 


Ai*ray for sinoothing 
0.5 of the filter ot\ier 
Wind flog 


f 

, 115 


1 


(46.72) 


i)lN 

rt*4 

JMP 

1*4 

IM 

1*4 

GAMMA 

0*8 


Array t» Ixj filtered 
latitudinal dimension 
Longitudinal dlmen^ on 
Filter parameters 


'^ItiGlNAL PAGE IS 
•>F P(X>R QUAMTV 






SUBROUTINE SIGTOP 


Argument list; A, B 


Description: Interpolates o-lovel quantities to the nandatory p levels; 

called by INCMN. 

Restrictions; Only those quantities designated by the NQ flag from ALTER2 
are Interpolated. 


Arguments and 
arrays passed 


Dimensions I S 


Description 


o-level profile 
p-level profile 


Gnninon areas 


DescrlDtion 


PHYSIC 

PRMANU 


Unstable profile count 
(not used) 

Data, flag from AIjTER 2 
Physical constants 
Pressure level constants 


External References 


Descriotion 


FIND 

PCAI 

FnjJN 

FINDHM 

HITEMP 

povmp 

RFLHUf4 

SBLI7.E 


Vertical interpolation (linear in log P) 
Competes o-level pressure 
Computes o-edge pressure 
Interpolates humidity profile 
Computes pressure level geopotential height 
Computes potential temperature 
Converts specific to relative humidity 
Stabilizes potential temperature profile 















SUBROUTINE SMOQC 


Argument list: PIN 


it 


I 


Description: Shuman filter applied to a horlaontal field; called ty l^dOOriH. 
Restrictions: None 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

P 


Description 

PIN 

R*4 

(46,72) 

D 

D 


Field to be filter 

Oonmcm areas 


length 



"^tescrlption 

OriRL 


060 



Model parameters 

Qttemal references 



Description 


1 



SHUM 


Shuman filter 






SUBlttJUTlNli SMOOTH 


Arguffient list. F, KS, TWINU 


Description, anoother for model fields after analysis called by 
INSEAL, INSURF and INUPAR. 

Smoottung flag KS siiould be 0 for no smoothing, 1 for Slavnan 1-diinetisional 
filter, 4 for 4-point shapiro filter, 8 for 8-point Shapiro filter^ or 9 for 
9-point amooth-desniooth . 

Restrictions: Set TWIND to .TRUE, for smoothing the wind field. 


Arguments and 
arrays passed 

Typo 

Dimensions 

S 

F 

A 

Description 

F 

R*4 

(46,72) 

X 

X 

X 

Field to be smootliod 

KS 

1*4 



X 


Snoothing flag 

TWIND 

L*4 



X 

X 

Wind flag 

Cbninon ai’eas 



Lengtii 


Description 

None 

External referetjces 



Description 

SMOOC 




Stiuman 

filter 

SHAi’n^ 




Shapiro filter 

$s^m) 




ST»ot h-desmootti 








SUBROUTINE SM2D 


Argument list: A, ISM 


Description: Five-point smoother not called. 

Restrictions: None 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

F 

A 

Description 

A 

R*4 

(46,72) 

X 

X 


Array to be smoothed 

ISM 

1*4 



X 


Number of applications 

Cbninon areas 




Description 

FAOCMf 


3300 



Work area 

CNTRL 


D60 



Model parameters 


SUBROUTINE SPCilUM 


Arguneat list 


H. T, P, L 


Description. Oanverts a relative humidity profile to a specific humidity 
profile, called by PTOSIG 


I 


I 

1 


I 

I 


Restrictions. None 


Argunents and 
arrays passed 

^ ■ 

Type 

Dimensions 

S 

1 

F 

A 

Description 

H 

R*4 

L 

X 

X 


Humidity profile 

T 


L 


X 

X 

Ibmperature profile 

P 

R*4 

L 


X 

X 

Pressure profile 

L 

1*4 



X 


Number of levels 

Ocninon areas 



Length 


Description 


None 


External references 

Description 

SATURN 

Saturation specific humidity 



SUBUnUTlNE SSCAN 


Argument list: JM, IM, RSCAN, CD, KNT 


Description: Station counting scanner; called by INSEAL, INSURE, and INUPAR 

This routine accumulates the number of enervations within a given radius 
about each grid point. 

Restrictions : None 


Arguments and 
arrays passed 

Type 

Dimensions 

S 

P 

A 

D^cription 

JM 

1*4 



X 


Latitudinal dimension 

IM 

1*4 



X 


Longitudinal dimension 

RSCAN 

R*4 



X 


Radius of scan in km 

CD 

R*4 

2 


X 


Earth coordinates of 
observations 

KNT 

1*4 

(JM,IM) 

X 

X 


Accumulator array 



NEXT Grid point scanner 

GRnXD Converts grid points to earth coordinates 

GCDIST Great circle distance calculation 



1 


r 

f 

► 

i- 

i 

} 

SUBROirriNE TAP24 

Description; Creates 
Restrictions: None 

Argument list: NTl, NT2 
LOG24 T^; called by ALTERS. 

1 

1 

Arguments and 

•Type 

DlmmislonB 

S 

D 

D 

DescrlptlcMi 

/* 

arrays passed 




■ 

■ 


t 

NTl 

1*4 





Input dl^ logical unit number 








(23) 


NT2 

1*4 


Oitput tape logical unit nunter 
(24) 

Odimcm 



tieacriptlon 

CNTRI. 


4DC 

Model parameters 

Model 


F2940 

Model arrays 

Auvr 


8 

Flag for LOG 24 tape 

External references 


DEHScrlptlon 


FVARIO 

GETEMP 

NEWTPE 


Provides D/A I/O 

Ootnputes inandator> level ten|)erature 
New output tape unit 24 






SUBIWGTINE TbaiP 


Argument list. JM, IM CD, F, PINT 


□oscripticMi. 
Restrictions . 

Porfonns a lioriksontal bilir ear interpolation. 
None 

Argunents and 
arrays passed 

Type 

Dijoenslons 

3 

F 

A 

Description 

JM 

IM 



X 


Latitudinal ditnension 

IM 

IM 



X 


Longitudinal dimension 

O) 

R*4 

2 


X 

X 

Earth coordinates 

K 

R*4 

(JM IM) 


X 


Grided field 

FIITT 

RM 


X 



Interpolated value 

Cbninon areas 



Length 


liescripti^ 

None 

External refei'ences 




Description 


aXJRlD 


Converts earth coordinates to grid coordinates 


SUBRCXJTINE TSPACT Argument list: P 


Description: Interpolates o-level temperatures to the surface. 


Restrictions: 

None 






Arguments and 
arrays passed 

TVpe 

Dimensi(XiB 

r"*“ ” 

S 

P 

A 

Description 

TSPACT 

P 

R*4 

R*4 

10 

X 

X 


Lapse rate returned 
o-level pressure 

Oomnon areas 


Length 



D^<:5*ipti^ 


Wiysicai constants 


PHYSIC 
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SUBIWTUTINE WINDPR Argunent list; NS 


Description: Uppor-alr wind check; called by ADPINS. *nils routine checks 
both conponents of the wind observations prior to the analysis. Erroneous 
winds are filled with 99999. 

Restrictions ; None 


Arguments and 
arrays passed 

Type 

Dimensions 

li- . 

S 



F 

A 

Description 

NS 

1*4 



X 


Number of reports 

Oomnon areas 



Length 


Descriptlcm 

(Same as INUPAR) 

External references 





Description 


LERI Single logical comparison 

CHKADP Lateral ched; 

CIJTADP Difference between observation and first guess check 

DIFFRS Difference between observation and first guess 


SUWWriNK WINiy2 


Arj^.»nt Ust; PR, PTC)P. T, TT 


Df?scrlption: Performs a linear in p Interpolation of a mandatory-lovei 

wind profile to tiie o levels, 

Restrl :tlons; Nont» 



r " - 


r — 

r - - * 

r •• *■ 


Am'uments and 

Type 

nirtKjnsloas 

s 

F 

A 

Description 

arrays paasjxl 




L - ^ 



PR 

R*4 



X 


Reference pressure 

PLOP 

R*4 



X 


Top pressun? 

T 

R*4 

12 


X 


Mandatory- level wind 

TT 

R*4 

9 

X 



0 -lovel wind 

Oonrion ai«as 



Length 


~ Description 

None 







External references 




Description 


Computes o-odRe pressures 


FIIJ.IN 


SUBROUTINE $SIKXX) 


ATiiUi«$nt i 1st QIN 


DuscrLption Nino-point anootU-dijsiaootti on a lilotjal two-diinuiisitmul 
cuiLud uy snuoni. 

tiustrictLons Mono 


Argumoiits and 1V|w DliffHisions ' 
arrays pasaud 


DuscrLption 


(46,72) 


Global model field 


(^nigm 

longth 

Description 

FAOXAI 

1U780) 

Work space 

CNTRL 

D50 

Forecast pfirametors 
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